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I.  INTRODUCTION 


The  aerodynamic  characteristics  of  standard  artillery  shell  from  subsonic 
to  supersonic  speeds  is  of  major  concern  in  the  design  of  new  shell  or  modifi¬ 
cations  to  existing  ones.  The  possibility  that  a  given  shape  may  have  to 
operate  throughout  a  range  of  Mach  numbers  requires  a  detailed  understanding 
of  the  flow  fields  associated  with  each  range.  Modern  computational  techni¬ 
ques  are  now  being  applied  to  projectile  shapes  and  the  ability  to  compute  the 
aerodynamics  of  shell  for  a  wide  range  of  Mach  numbers  and  Reynolds  numbers  is 
becoming  a  reality.  Significant  accomplishments  have  been  made  in  the  super¬ 
sonic  regime1:  static  and  Magnus  force  coefficients  have  been  computed  for 
standard  projectile  configurations  and  experimental  data  is  generally  avail¬ 
able  which  shows  good  comparisons  with  theory.  Transonic  flow,  however, 
presents  a  new  complexity  for  computational  analysis.  The  formation  of  shock 
waves,  imbedded  in  the  flow  field  near  surface  discontinuities,  produces  a 
severe  change  of  the  aerodynamic  coefficients  such  as  drag  and  pitching 
moments.  For  example,  the  drag  of  a  projectile  shape  has  been  found  to  change 
by  as  much  as  100%  through  a  Mach  number  range  of  .95  to  .97.  A  change  of 
this  magnitude  in  the  aerodynamics  makes  it  essential  to  understand  and 
compute  the  features  of  the  flow  field  which  contribute  to  this  effect. 

The  transonic  flow  field  around  a  projectile  presents  a  difficult  and 
interesting  problem;  the  general  features  of  the  flow  pattern  are  illustrated 
in  Figure  1,  which  is  a  shadowgraph  of  a  typical  projectile  flying  at 
M.  =  .946.  The  two  surface  shape  discontinuities  create  local  expansion  zones 
terminated  by  shock  waves  through  which  the  boundary  layer  must  pass,  thereby 
creating  a  pair  of  vi scous-i nvi scid  interaction  regions.  The  existance  of 
these  regions  together  with  the  cylindrical  geometry  of  projectiles  clearly 
produces  a  composite  i nvi sci d- vi scous  flow  field  problem  not  encountered  in 
the  design  of  normal  winged  aircraft. 

Although  some  comprehensive  transonic  flow  field  investigations  have  been 
made  on  transonic  wings2’3,  comparable  efforts  have  not  been  made  for  bodies 
of  revolution  particularly  if  spinning.  Existing  studies  have  been  primarily 
concerned  with  surface  pressure  distributions  and  inviscid  flow  predictions 
with  only  sparse  consideration  of  the  boundary  layer  and  transonic  viscous- 


1.  Sturek,  W.  B.,  et  at.,  "Computations  of  Magnus  Effects  for  a  Yawed , 
Spinning  Body  of  Revolution  }'*  AIM  Journal,  Vol .  16,  No.  7,  July 
1978,  pp.  687-692. 

2.  Rose,  W.  C . ,  and  Segniner,  A.,  "Calculation  of  Transonic  Flow  Over 
Supercritical  Airfoil  Sections,"  AIAA  Journal  15,  August  1978 ,  pp.  514- 
519. 

3.  Sobieczky,  M.,  and  Stanewsky ,  E.,  "The  Design  of  Transonic  Airfoils  Under 
Consideration  of  Shock  Wave-Boundary  Layer  Interaction, "10  ICAS  Congress 
Paper,  October  1976  (see  also  E.  Stanewsky  and  H.  Zimmer,  Zeit.  fur 
Flugwissenschaften  23,  Heft  7/8,  1975). 
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inviscid  interaction  effects  throughout  the  flow  field4’5.  Consequently,  a 
comprehensive  fundamental  investigation  of  the  combined  i nvi sci d-vi scous  flow 
field  problem  on  transonic  artillery  projectiles  was  undertaken  at  BRL  with 
the  objective  of  developing  an  appropriate  composite  theoretical  treatment 
including  the  shock-boundary  layer  interaction  effects,  together  with  validat¬ 
ing  comparisons  with  experiment.  The  present  report  describes  the  results  of 
this  study  for  the  case  of  projectile  flight  at  zero  angle  of  attack,  with 
emphasis  on  a  description  of  the  theoretical  work  involved  in  treating  the 
interaction  problem.  Since  the  approach  taken  has  been  to  extend  Inger's 
theory6  of  non-separating  transonic  shock-turbulent  boundary  layer  interac¬ 
tions  to  spinning  bodies  of  revolution,  we  first  highlight  the  2-D  version  of 
this  theory  in  Section  2.  The  desired  extension  of  the  theory  is  then  given 
in  Section  3,  followed  in  Section  4  by  a  description  of  how  the  local  interac¬ 
tion  theory  is  combined  with  inviscid  transonic  and  turbulent  boundary  layer 
prediction  methods  to  form  a  composite  global  i nvi sci d-vi scous  flow  field  pre¬ 
diction  code  for  the  entire  body.  Finally,  Section  5  presents  comparisons  of 
this  theory  with  both  experimental  data  and  numerical  solutions  based  on  the 
thin-layer  Navier-Stokes  equations. 


II.  OUTLINE  OF  THE  LOCAL  INTERACTION  THEORY  FOR  TWO-DIMENSIONAL  FLOWS 

A.  Rationale  of  Non-Asymptotic  Triple  Deck  Approach 

It  is  well-known7  that  when  separation  occurs,  the  disturbance  flow  pat¬ 
tern  associated  with  a  nearly-normal  shock-boundary  layer  interaction  is  a 
very  complicated  one  involving  a  bifurcated  shock  pattern,  whereas  the  unsepa¬ 
rated  case  pertaining  to  turbulent  boundary  layers  up  to  roughly  Mi  «  1.3  has 

instead  a  much  simpler  type  of  interaction  pattern  which  is  more  amenable  to 
analytical  treatment  (see  Figure  2).  With  some  judicious  simplifications,  it 
is  possible  to  construct  a  fundamentally-based  approximate  theory  of  the 


4.  Wu,  J.  C.,  Moulden,  T.H.,  and  Uchiyama,  N.,  "Aerodynamic  Performance  of 
Missile  Configurations  at  Transonic  Speeds  Including  the  Effects  of  a  Jet 
Plume , "  U.S.  Army  Missile  Command  Technical  Report  RD-76-23,  Redstone 
Arsenal ,  Alabama ,  March  1976. 

5.  Reklis ,  R.  P.,  Sturek ,  W.  B.,  and  Bailey ,  F.  R.,  "Computation  of  Transonic 
Flow  Past  Projectiles  at  Angle  of  Attack ,"  AIAA  11th  Fluid  and  Plasma 
Dynamics  Conference ,  AIAA  Paper  No.  78-1182,  July  1978. 

6.  Inger,  G.  R.,  "Upstream  Influence  and  Skin  Friction  in  Non- Separating 
Shock  Turbulent  Boundary  Layer  Interactions ,  "  AIAA  Pamper  80^1411^ 

Snowmass,  Colorado,  July  1980. 

7.  Ackeret,  J.,  Feldman,  F.,  and  Rott,  N.,  "Investigations  of  Compression 
Shocks  and  Boundary  Layers  in  Gases  Moving  at  High  Speed,"  NACA  TM-1113, 
January  1947. 
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problem,  as  documented  in  detail  in  Reference  6.  For  purposes  of  orientation 
and  completeness,  a  brief  outline  of  this  theory  will  be  given  here. 

We  consider  small  disturbances  of  an  arbitrary  incoming  turbulent  boundary 
layer  due  to  a  weak  external  shock  and  examine  the  detailed  perturbation  field 
within  the  layer.  We  purposely  employ  a  non-asymptotic  triple-deck  flow  model 
patterned  in  some  ways  after  Lighthill's  approach8  because  of  its  essential 
soundness  and  adaptability  to  further  improvement,  because  of  its  similarity 
to  related  types  of  multiple-deck  approaches  that  have  proven  highly  success¬ 
ful  in  treating  turbulent  boundary  layer  response  to  strong  known  adverse 
pressure  gradients9,  and  because  of  the  large  body  of  turbulent  boundary  layer 
interaction  data  that  supports  the  predicted  results  in  a  variety  of  specific 
problems6.  At  high  Reynolds  numbers  it  has  been  established10’* 11’12  that  the 
local  interaction  disturbance  field  in  the  neighborhood  of  the  impinging  shock 
organizes  itself  into  three  basic  layered-regions  or  "decks"  (Figure  3):  1) 

an  outer  region  of  potential  inviscid  flow  above  the  boundary  layer,  which 
contains  the  incident  shock  and  interactive  wave  systems;  2)  an  intermediate 
deck  of  frozen  shear  stress-rotational  inviscid  disturbance  flow  occupying  the 
outer  90%  or  more  of  the  incoming  boundary  layer  thickness;  3)  an  inner 
shear-disturbance  sublayer  adjacent  to  the  wall  which  accounts  for  the  inter¬ 
active  skin  friction  perturbations  (and  hence  any  possible  incipient  separa¬ 
tion)  plus  most  of  the  upstream  influence  of  the  interaction.  The  "forcing 
function"  of  the  problem  here  is  thus  impressed  by  the  outer  deck  upon  the 
boundary  layer;  the  middle  deck  couples  this  to  the  response  of  the  inner  deck 
but  in  so  doing  can  itself  modify  the  disturbance  field  to  some  extent,  while 
the  slow  viscous  flow  in  the  thin  inner  deck  reacts  very  strongly  to  the 
pressure  gradient  disturbances  imposed  by  these  overlying  decks.  This  triple 


8.  Lighthill ,  M.  J .,  "On  Boundary  Layers  and  Upstream  Influence;  II. 
Supersonic  Flow  Without  Separation ,"  Pmcr  Rnual  8oc.  A  21  7  r  1953 ,  pp. 
478-507. 

9.  Stratford,  B.  S.,  "The  Prediction  of  Separation  of  the  Turbulent  Boundary 
Layer,"  Jqkt*.  Eluid  MMh*  5,  pp.  1-16 ,  1959. 

10.  Inger,  G.  R.,  and  Mason,  W.  H .,  "Analytical  Theory  of  Transonic  Normal 
Shock-Boundary  Layer  Interaction,"  AIM  Eap£2  75-831  June  1975 
(abbreviated  version  published  in  AIM.  Iqujl.  1£>  pp*  1266-72,  September 
1976). 

11.  Melnik,  R.  E.,  and  Grossman,  B.,  "Analysis  of  the  Interaction  of  a  Weak 
Normal  Shock  Wave  with  a  Turbulent  Boundary  Layer, "  AIM  E2££2l  74-59 H . 
June  1974. 

12.  Adamson,  T.  C.,  and  Feo,  A.,  "Interaction  Between  a  Shock  Wave  and  a 
Turbulent  Layer  in  Transonic  Flow,"  SIAM  J£U£m1jl  Appl.  Math  29.  July 
1975,  pp.  121-144. 
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deck  structure  also  has  been  employed  in  the  theoretical  studies  of  Gadd13, 
Lighthill8,  Stratford9,  Honda14  and  others,  and  has  been  verified  by  a  large 
body  of  experimental  evidence  and  recent  numerical  studies  with  the  full 
Navier-Stokes  equations15.  The  essential  correctness  of  this  model  is  further 
supported  by  its  success  in  related  boundary  layer  perturbation  problems 
involving  vi scous-i nvi scid  interactions,  turbulent  boundary  layer  response  to 
sudden  changes  in  surface  roughness  or  pressure  gradient,  and  flow  past 
various  kinds  of  surface  distortions  including  skin  friction  measuring  devices 
(see  Reference  6) . 

While  there  is  general  agreement  about  the  validity  of  the  triple  deck 
approach  over  a  wide  range  of  Reynolds  number  and  the  well-known  qualitative 
differences  in  interactive  response  between  laminar  and  turbulent  flow  (e.g., 
the  much  smaller  upstream  influence  and  larger  separation-resistance  of  the 
later),  questions  have  been  raised  concerning  (a)  the  relative  interactive 
importance  of  the  inner  shear-disturbance  deck  and  (b)  the  accuracy  of  delib¬ 
erately  using  a  non-asymptoti c  treatment  of  the  details  within  the  boundary 
layer.  Regarding  the  first,  we  note  that  while  asymptotic  (Re^  +  °°)  theory 

predicts  an  exponenti al ly-smal 1  thickness  and  displacement  effect  contribution 
of  the  inner  deck9’10,  this  is  not  apparently  true  at  ordinary  Reynolds 
numbers,  where  many  analytic  and  experimental  studies  have  firmly  established 
that  this  deck,  although  indeed  very  thin,  still  contributes  significantly  to 
the  overlying  interaction  and  its  displacement  thickness  growth6.  Bolstered 
by  these  facts,  plus  the  unanimous  conclusion  reached  in  the  detailed  reviews 
by  Green16,  Rose,  Murphy  and  Watson17,  and  Hankey  and  Holden18  that  the 
viscous  sublayer  is  an  important  component  of  turbulent  interaction  problems, 
we  take  the  point  of  view  here  that  the  inner  deck  is  in  fact  significant  at 


13.  Gadd ,  G.  E.,  "Interactions  Between  Wholly  Laminar  or  Wholly  Turbulent 
Boundary  Layers  and  Shock  Wave  Strong  Enough  to  Cause  Separation ,  "  Jour . 
of  the  Aeronaut.  Sci.  20,  November  1953,  p.  729. 

14.  Honda,  M.,  "A  Theoretical  Investigation  of  the  Interaction  Between  Shock 
Waves  and  Boundary  Layers, "  Jour.  Aero/Space  Sci.  25,  November  1958,  pp. 
667-677. 

15.  Hankey,  W.,  and  Shang,  J .,  'Numerical  Solution  of  the  Navier-Stokes 
Equations  for  Supersonic  Turbulent  Flow  over  a  Compression  Ramp, "  AIAA 
Paper  75-3,  Pasadena,  Jan.  1975. 


16.  Green,  J.  E.,  "Interactions  Between  Shock  Waves  and  Boundary  Layers,"  in 
Progress  of  Aero.  Sci.,  Vol.  11,  Pergamon,  N.Y.,  1965,  pp.  319. 

17.  Rose,  W.  C.,  Murphy,  J.  D .,  and  Watson,  E.  C.,  "Interaction  of  an  Oblique 
Shock.  Wave  with  a  Turbulent  Boundary  Layer,  "  AIAA  Jour.  7,  December  1969, 

pp.  2211-2221 . 

18.  Hankey,  W.  L.,  and  Holden,  M.  S.,  "Two-Dimensional  Shock  Wave-Boundary 
Layer  Interactions  in  High  Speed  Flow,  AGARDograph  203,  June  1974. 
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the  Reynolds  numbers  of  practical  interest.  In  this  regard,  we  re-emphasize 
that  this  deck  contains  all  of  the  skin  friction  and  incipient  separation 
effects  in  the  interaction,  which  alone  are  suficient  reasons  to  examine  it  in 
detail.  Regarding  (b),  it  is  pointed  out  that  application  of  Re  ■*■  ®  asymp- 

totic  theory  results  (no  matter  how  rigorous  in  this  limit)  to  ordinary 

Reynolds  numbers  is  itself  an  approximation  which  may  be  no  more  accurate, 
(indeed  perhaps  less  so)  than  a  physically  well-constructed  non-asymptoti c 
theory.  Direct  extrapolated-asymptotic  versus  non-asymptoti c  theory  compari¬ 
son  have  definitely  shown  this  to  be  the  case  for  laminar  flows  (especially  as 
regards  the  skin  friction  aspect19)  and  the  situation  has  been  shown  to  be 
possibly  even  worse  in  turbulent  flow6.  For  example,  the  asymptotic  first 

order  theory  formally  excludes  both  the  streamwise  interactive  pressure 

gradient  effect  on  the  shear-disturbance  deck  and  both  the  normal  pressure 
gradient  and  so-called  "streamline  divergence"  effects  on  the  middle  deck; 

however,  physical  considerations  plus  experimental  observations  and  recent 
comparative  numerical  studies20’21  suggest  that  these  effects  may  in  fact  be 
significant  at  practical  Reynolds  numbers  and  should  not  be  neglected.  Of 
course,  second  order  asymptotic  corrections  can  be  devised  to  redress  this 
difficulty  but,  as  Neyfeh  and  Regab22  have  shown,  run  the  risk  of  breaking 
down  even  worse  when  extrapolated  to  ordinary  Reynolds  numbers.  In  the 
present  work,  we  avoid  these  problems  by  using  a  deliberately  non-asymptoti c 
triple-deck  model  appropriate  to  realistic  Reynolds  numbers  that  includes  the 
inner  deck  pressure  gradient  terms  plus  the  middle  deck  3p/3y  and  streamline 
divergence  effects,  along  with  some  simplifying  approximations  that  render  the 
resulting  theory  tractible  from  an  engineering  standpoint.  With  this  view¬ 
point  in  mind,  we  now  examine  in  more  detail  the  nature  of  the  disturbance 
flow  problem  in  each  of  the  three  basic  decks. 

B,  Formulation  of  the  Disturbance  Problem  in  Each  Deck 


1.  Outer  Potential  Flow  Region.  Assuming  that  the  incident  shock  and  its 
reflection  system  are  weak  with  isentropic  non-hypersonic  flow,  we  have  here  a 


19.  Burggraf ',  0.  R .,  "Asymptotic  Theory  of  Separation  and  Attachment  of  a 
Laminar  Boundary  Layer  on  a  Compression  Ramp , "  in  AGARD  CP-168  -  Flow 
Separation ,  1975. 

20.  Werle ,  M.,  and  Berthe,  S.  D.,  "Application  of  an  Interacting  Boundary 
Layer  Model  to  the  Supersonic  Turbulent  Separation  Problem, "  University 
of  Cincinnati  Report  AFL  76-4-21,  August  1976. 

21.  Le  Balleur,  J.  C.,  Peyret,  R.,  and  Vivand,  H.,  "Numerical  Studies  in  High 
Reynolds  Humber  Aerodynamics , "  in  Computers  and  Fluids,  Vol.  8,  Pergamon, 
1980,  pp.  1-30. 

22.  Ragab,  S.  A.,  and  Hay f eh,  A.  H.,  "A  Second  Order  Asymptotic  Solution  for 

Laminar  Separation  of  Supersonic  Flows  Past  Compression  Ramps, "  AIAA  78- 
1132,  1978.  - 
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small  disturbance  potential  inviscid  motion  imposed  upon  the  undisturbed 

uniform  flow  U.  outside  the  boundary  layer: 
ue 


[M  2  -1  +  (y  +  1)  u' 


1  3u 1 
J  3x 


3  V 1 

3y 


(i) 


3v'/9x  =  3  U  1  /  3y 


(2) 


2_  i 


3  P 

3x2 


0 


(3) 


where  the  third  term  within  the  square  brackets  of  Eqs.  (1)  and  (3)  is  signif¬ 
icant  in  the  transonic  regime  1  <  MQ  <1.1  and  automatically  includes  the 

6 

supersonic-subsonic  jump  conditions  to  this  order  of  approximation23.  Since  a 
variety  of  efficient  analytical  or  numerical  methods  are  presently  available 
to  solve  this  system  in  either  transonic  flow  or  in  purely  supersonic  flow  (in 
which  case  Eqs.  1-3  further  reduce  to  an  Ackert-type  problem),  we  assume  that 
such  a  solution  may  be  carried  out  for  all  x  on  the  upper  region  y  >  <SQ  sub¬ 
ject  to  the  usual  far-field  conditions  as  y  ».  The  remaining  disturbance 
boundary  condition  that  must  be  supplied  along  y  =  <$Q  then  couples  this  solu¬ 
tion  to  the  underlying  double-deck:  it  requires  that  the  outer  disturbance 
flow  pertain  to  an  effective  streamline  shape  (relative  to  the  wall)  defined 
by  the  total  interactive  displacement  effect  of  the  inner  decks.  To  insure 
physi cal ly-smooth  matching  along  this  outer-inner  interface,  then,  we  require 
both  v'/U  and  p'  to  be  continuous  with  their  middle  deck  counterparts  along 
°e 


2.  Middle  Rotational -Disturbance  Flow  Deck.  This  layer  contributes  to 
and  transmits  the  di splacement  effect,  contai ns  the  boundary  layer  lateral 
pressure  gradient  due  to  the  interaction  and  carries  the  significant  influence 
of  the  incoming  boundary  layer  profile  shape.  Our  analysis  of  this  layer 
rests  on  the  key  simplifying  assumption  that  for  non-separating  interactions 
the  turbulent  Reynolds  shear  stress  changes  are  small  and  have  a  negligible 


23.  Murrmn ,  E.  M.,  and  Cole,  J.  D .,  "Calculation  of  Plane  Steady  Transonic 
Flow,  "  AIM  Jour.  9,  January  1971,  pp.  114-121. 
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back  effect  on  the  mean  flow  properties  along  the  interaction  zone;  hence  this 
stress  can  be  taken  to  be  "frozen"  along  each  streamline  at  its  appropriate 

value  in  the  undisturbed  incoming  boundary  layer.  This  approximati on ,  like¬ 
wise  adopted  by  a  number  of  earlier  investigators  with  good  results,  is  sup¬ 
ported  not  only  bv  asymptotic  analysis24  but  especially  by  the  results  of 

Rose's  several17’2^  detailed  experimental  studies  of  a  non-separati ng  shock- 
turbulent  boundary  layer  interaction  which  showed  that  over  the  short-ranged 
interaction  length  straddling  the  shock  the  pressure  gradient  and  inertial 
forces  outside  a  thin  layer  near  the  wall  are  at  least  an  order  of  magnitude 

larger  than  the  corresponding  changes  in  Reynolds  stress.  Furthermore,  there 

is  a  substantial  body  of  related  experimental  results  on  turbulent  boundary 
layer  response  to  various  kinds  of  sudden  perturbations  and  rapid  pressure 
gradients  which  also  strongly  support  this  view  (see  Ref.  6).  These  studies 
unanimously  confirm  that,  at  least  for  non-separati  ng  flows,  significant  local 
Reynolds  shear  stress  disturbances  are  essentially  confined  to  a  thin  sublayer 
within  the  Law  of  the  Wall  region  (see  below)  where  the  turbulence  rapidly 
adjusts  to  the  local  pressure  gradient,  while  outside  in  the  Law  of  the  Wake 
region  the  turbulent  stresses  respond  very  slowly  and  remain  nearly  frozen  at 
their  initial  values  far  out  of  local  equilibrium  with  the  wall  stress. 

Confining  attention,  then,  to  the  short  range  local  shock  interaction  zone 
where  the  aforementioned  "frozen  turbulence"  approximation  is  applicable,  the 
disturbance  field  caused  by  a  weak  shock  is  one  of  small  rotational  inviscid 
perturbation  of  the  incoming  non-uniform  turbulent  boundary  layer  profile 
M0(y),  governed  by  the  equations 


^  , /  x  1-M  ^(y)  3(p ' /p  ) 

3  .  r  v  (x,y)  i  o  w/  /hoy 

*  l^r]  =  7T7^  • 


(4) 


3u 1  3p 1  /  3x  _  ^o  #  vl 

"  P0(y)  UQ(y)  dy  *  UQ 


24.  Yajnik,  A "Asymptotic  Theory  of  Turbulent  Wall  Boundary  Layer  Flows," 
JFM  42j  1970 ,  pp.  411-427. 

25.  Rose ,  W.  C .,  and  Childs ,  M.  E .,  "Reynolds  Shear  Stress  Measurements  in  a 
Compressible  Boundary  Layer  within  a  Shock  Wave-Induced  Adverse  Pressure 
Gradient ,  ”  JFM  65,  1,  1974 ,  pp.  177-188. 
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(6) 


9  dM„  a. 

1 _ 2  i£_  + 

MQ  dy  3y 


1-M. 


2  u'M 
U 


=  0 


where  Eq.  (4)  is  a  result  of  the  combined  parti cle-i sentropi c  continuity,  x- 
momentum  and  energy  equations.  It  is  noted  that,  consistent  with  the  assumed 
short  range  character  of  the  interaction,  the  streamwise  variation  of  the 
undisturbed  turbulent  boundary  layer  properties  that  would  occur  over  this 
range  are  neglected,  taking  U0(y),p0(y)  and  M0(y)  to  be  arbitrary  functions  of 

y  only  with  6  ,  6  *  and  x  as  constants.  Now  Eq.  (6)  is  a  generalization  of 
0  0  w  _ 

0 

Lighthill's  well-known  pressure  perturbation  equation  for  non-uniform  flows8 
which  includes  a  non-linear  correction  term  for  possible  transonic  effects 
within  the  boundary  layer  including  the  diffracted  impinging  shock  above  the 
sonic  level  of  the  incoming  boundary  layer  profile.  Excluding  the  hypersonic 
regime,  Eqs.  (4)-(6)  therefore  apply  to  a  wide  range  of  initially  supersonic 
external  flow  conditions  and  the  complete  speed  range  across  the  boundary 
layer  except  at  the  singular  point  MQ  -*■  0  (which  we  avoid  by  consideration  of 

the  inner  deck  as  shown  below).  In  particular,  use  of  Eq.  (6)  provides  an 
account  of  any  lateral  pressure  gradient  that  develops  across  the  interacting 
boundary  layer. 


As  is  the  case  with  the  outer  deck,  a  variety  of  analytical  or  numerical 
methods  may  be  used  to  solve  this  middle  deck  disturbance  problem  (see  e.g. , 
Refs.  8,  10,  17).  Whatever  the  method  chosen,  we  imagine  that  it  provides  at 
each  streamwise  station  x  an  evaluation  of  the  disturbance  pressure  distribu¬ 
tion  p'(x,y);  then  .y-integration  of  Eq.  (4)  gives 


v‘(*»y)  = 
V7T 


(x,y 


)  +  _ 

weff  8x 


=  0 


1-M02(y) 

M02(y) 


(7) 


where  y  >  0  is  the  effective  wall  shift  or  displacement  height  associated 
weff 

with  the  inner  deck  defined  such  that  the  inviscid  v1  (x,yw  1  and  hence 

"eff 

9p'/3y(x,y  )  both  vanish  (see  below).  Equation  (7)  provides  the  disturb- 

weff 

ance  streamline  slope  distribution  across  the  boundary  layer  at  any  streamwise 
station,  and  its  value  at  y  =  6q  yields  the  total  streamline  displacement 

effect  of  the  two  inner  decks  (the  lower  limit  on  the  integral  being  the  inner 
deck  contribution).  We  then  may  obtain  the  corresponding  total  displacement 
thickness  growth  along  the  interaction  by  streamwi se-quadrature  of  the  pertur¬ 
bation  boundary  layer  continuity  equation  integral;  this  yields  to  first 
perturbation 
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A<5*( x)  “  /  ° 


<5 


(8) 


eff 


3.  The  Inner  Shear-Disturbance  Layer.  This  very  thin  inner  deck  contains 
the  si gni fi cant  vi scous  ancl  turbulent  shear  stress  disturbances  due  to  the 
interaction,  plus  the  small  upstream  influence  and  an  important  contribution 
to  the  viscous  displacement  effect.  In  fact  it  lies  well  within  the  Law  of 
the  Wall  region  of  the  incoming  turbulent  boundary  layer  profile  and  also  (for 
the  high  Reynolds  numbers  of  interest  here)  below  the  sonic  level  of  the  pro¬ 
file  as  our  resulting  theory  indeed  confirms  aposteriori.  The  original  work 
of  Lighthill8  and  others  treated  the  problem  by  further  neglecting  the  turbu¬ 
lent  stresses  altogether  and  considering  only  the  laminar  sublayer  effect; 
while  this  greatly  simplifies  the  problem  and  yields  an  elegant  analytical 
solution,  the  results  can  be  significantly  in  error  at  high  Reynolds  numbers 
and  cannot  explain  (and  indeed  conflicts  with)  the  ultimate  asymptotic 
behavior  pertaining  to  the  Re£  +  «  limit.  The  present  theory  remedies  this  by 

extending  Lighthill 's  approach  to  include  the  entire  Law  of  the  Wall  region 
turbulent  stress-effects;  the  resulting  general  shear-disturbance  sublayer 
theory  provides  a  non-asymptoti c  treatment  which  encompasses  the  complete 
range  of  Reynolds  numbers.  It  is  noted  in  this  connection  that  our  considera¬ 
tion  of  the  entire  Law  of  the  Wall  combined  with  the  use  of  the  effective 
inviscid  wall  concept  to  treat  the  inner  deck  displacement  effect  eliminates 
the  need  for  the  “blending  layer"11  that  is  otherwise  required  to  match  the 
disturbance  field  in  the  laminar  sublayer  region  with  the  middle  inviscid 
deck;  except  for  higher  order  derivative  aspects  of  asymptotic  matching,  our 
inner  solution  effectively  includes  this  blending  function  since  it  imposes  a 
boundary  condition  of  vanishing  total  shear  disturbance  at  the  outer  edge  of 
the  deck.  In  addition,  our  retention  of  the  explicit  disturbance  pressure 
gradient  term  for  the  inner  deck  not  only  provides  the  correct  physics  at 
practical  Reynolds  numbers  but  also  correctly  models  the  situation  near  sepa¬ 
ration  (t  >0)  where  this  term  becomes  of  dominant  importance. 

To  facilitate  a  tractible  theory,  we  retain  only  the  main  physical  effects 
by  introducing  the  following  simplifying  assumptions,  (a)  The  incoming  bound¬ 
ary  layer  is  free  from  any  post-transi ti onal  memory  or  low  Reynolds  number 
effects  and  its  Law  of  the  Wall  region  is  characterized  by  a  constant  total 
(laminar  plus  turbulent  eddv)  shear  stress  and  a  Van  Driest-Cebeci  type  of 
damped  eddy  viscosity  model  ^  .  This  model  is  known  to  be  a  good  one  for  a 
wide  range  of  upstream  non-separating  boundary  layer  flow  histories,  (b)  For 
the  weak  incident  shock  strengths  of  present  interest,  the  sublayer  disturb¬ 
ance  flow  is  assumed  to  be  a  small  perturbation  upon  the  incoming  boundary 


26.  Cebeci ,  T.,  and  Bradshaw,  P.,  ''Mormrvbmn  Transfer  in  BoundaanL  Lavers. 
McGraw-Hill/Hemisphere,"  Wash.,  D.  C.,  1977,  p.  365. 
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layer;  in  the  resulting  linearized  disturbance  equations,  however,  al 1  the 
physi cal ly-important  effects  of  streamwise  pressure  gradient,  streamwise  and 
vertical  acceleration,  and  both  laminar  and  turbulent  disturbance  stresses  are 
retained.  Although  the  resulting  theory  necessarily  becomes  inaccurate  near 
separation,  it  provides  a  valuable  physical  insight  to  the  interactive  physics 
close  to  the  wall  for  non- separating  flow  and  a  firm  basis  for  subsequent 
improvement.  Moreover,  it  can  be  shown  that  the  form  of  the  particular  set  of 
linear  equations  used  here  is  in  fact  unaltered  by  non-linear  effects;  hence 
even  the  quantitative  accuracy  of  the  present  theory  is  expected  to  be  good 
until  rather  close  to  separation,  (c)  For  adiabatic  flows  at  1 ow-to-moderate 
external  Mach  numbers,  the  undisturbed  and  perturbation  flow  Mach  numbers  are 
both  quite  small  within  the  shear  disturbance  sublayer;  consequently  the 
treatment  of  compressibility  effects  therein  can  be  greatly  simplified  without 
any  significant  error.  Thus,  the  influence  of  the  density  perturbations  on 
the  sublayer  disturbance  flow  may  be  neglected  altogether,  while  the  corre¬ 
sponding  modest  compressibility  effect  on  the  Law  of  the  Wall  portion  of  the 
undisturbed  profile  is  quite  adequately  treated  by  the  Eckert  reference 
temperature  method27  wherein  incompressible  relations  are  used  based  on  wall 
recovery  temperature  properties.  Excluding  hypersonic  flow  this  is  equivalent 
in  accuracy  to  (but  easier  than)  the  use  of  Van  Driest's  compressible  Law  of 
the  Wall  profile28,  (d)  The  turbulent  fluctuations  and  the  small  interactive 
disturbances  are  assumed  uncorrelated  in  both  the  lower  and  middle  decks,  (e) 
The  thinness  of  the  inner  deck  allows  the  boundary  layer-type  approximation  of 
neglecting  its  lateral  pressure  gradient;  the  wall  pressure  distribution  p^(x) 

is  taken  equal  to  the  overlaying  pressure  perturbation  field  along  the  bottom 
of  the  middle  deck. 

Under  these  assumptions,  the  disturbance  field  is  governed  by  the  follow¬ 
ing  continuity  and  momentum  equations: 


(9) 


(10) 


(11) 


3p'/9y  =  o  ;  p‘  =  p'(x)  “  p^(x) 


27 .  Burggraf,  0 .  R.,  "The  Compressibility  Transformation  and  the  Turbulent 
Boundary  Layer  Equation ,  "  Jour,  of  the  Aerospace  Sci.  29,  1962 ,  pp.  434- 
439. 

28.  Van  Driest ,  E.  R.,  "Turbulent  Boundary  Layers  in  Compressible  Fluid," 
Jour,  of  the  Aergr]gut_.  Sci.  18 .  March  1951,  pp.  145-160. 
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where  p  and  v  are  evaluated  at  the  adiabatic  wall  recovery  temperature  and 
w  w 

0  0 

where  it  should  be  noted  that  the  kinematic  eddy  viscosity  perturbation  Sj'  is 

being  taken  into  account.  The  corresponding  undisturbed  turbulent  boundary 
layer  Law  of  the  Wall  profile  UQ(y)  is  governed  by 


x  (y)  =  const, 
o 


u 


w 


dU 

P  eT  (y)3  ~a — ~ 

WQ  T0vy;j  dy 


(12) 


where  according  to  the 
y+  =  (y  /xTTTD/v. 


Van  Driest-Cebeci  eddy  viscosity  model  with 


eT  =  [.41y(l-e-y+/A)]2-|^  (13A) 

which  yields  for  non-separating  flow  disturbances  that 

eT  =  [.41y(l-e'y  /A)]2-xr  (13B) 

o  y 


,  3u'/3y 
ldUQ/dy 


)e 


T 


o 


(13C) 


Here,  A  is  the  so-called  Van  Driest  damping  "constant";  we  used  the  commonly- 
accepted  value  A  =  26  although  it  is  understood  that  a  larger  value  may 
improve  the  experimental  agreement  in  regions  of  shock-boundary  layer  interac¬ 
tion.29  Substituting  (13c)  into  (10)  we  thus  have  the  disturbance  momentum 
equati on 


3u 1 
3x 


+  v 


dUo 

w 


(Pw  _1) 


3p 1 

9X 


■jy 


+  it 


)  3U  >1 
T  >  3y  J 
o  J 


(14) 


from  which  we  see  that  inclusion  of  the  eddy  viscosity  perturbation  has 
exactly  doubled  the  turbulent  shear  stress  disturbance  term. 


29.  Jobe ,  C.  E.,  and  Harikey ,  W.  L.,  "Turbulent  Boundary  Layer  Calculations  in 
Adverse  Pressure  Gradient  Flows , "  AIAA  Pager  80-0136,  Pasadena,  January 
1980. 
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We  seek  to  solve  Eqs.  (9)  and  (14)  subject  to  the  wall  boundary  conditions 
Uo(0)  =  u'(x,o)  =  v'(x,o)  =  0  plus  an  initial  condition  u'(-°°,y)  =  0  requiring 

that  all  interactive  disturbances  vanish  far  upstream  of  the  impinging  shock. 
Furthermore,  at  some  distance  6^^  sufficiently  far  from  the  wall,  u1  must  pass 

over  to  the  inviscid  solution  ujnv  along  the  bottom  of  the  middle  deck,  this 

latter  being  governed  by  Eq.  (9)  plus 


9u . 
inv 

3x 


+  v! 


l  nv 


dU 

-r-2-  +  (p 

dy  v  w 


3p 
_ v 

3x 


«  0 


(15) 


with  defined  as  the  height  where  the  total  shear  disturbance  (proportional 
to  3u'/3y)  of  the  inner  solution  vanishes  to  a  desired  accuracy. 

Following  Lighthill8,  it  proves  convenient  to  convert  the  foregoing 
problem  into  one  involving  the  normal  disturbance  velocity  field  v'(x,y). 
Differentiating  Eq.  (14)  w.r.t.  x,  substituting  Eq.  (9)  so  as  to  eliminate  u' 
and  then  differentiating  the  result  w.r.t.y  so  as  to  eliminate  p^  by  virtue  of 

Eq.  (11),  one  thus  obtains  the  following  fourth-order  equation  for  v ' : 


J 

3x 


"  dy2 


32V 1  i 

3y^~ 


(16) 


This  differential  equation  contains  a  three-fold  influence  of  the  turbulent 
flow:  the  profile  U0(y) ,  its  curvature  d^/dy2  (non-zero  outside  the  laminar 

sublayer)  and  a  new  eddy  disturbance  stress  term  2 e-j.  .  Equation  (16)  is  to  be 

o 

solved  together  with  (11)  and  (13b)  subject  to  the  wall  boundary  conditions 
v'(x,o)  =  3v'/3y(x,o)  =0.  A  third  condition  involving  v'  is  obtained  by 
satisfying  the  x-momentum  equation  (14)  at  the  wall;  when  this  is  differenti¬ 
ated  w.r.t.  x  and  Eq  (9)  used  together  with  the  fact  that  e_  ■>■  0  at  y  =  0  we 
obtain  the  non- homogeneous  condition 


33  v1 

3y3 


(x,o) 


d2p 1 
Hw 

dx2 


(17) 


The  fourth  boundary  condition  is  the  v1  equivalent  of  the  outer  inviscid 
matching  requirement  (15),  which  yields  v'<x-5SL>  '  ’'inv<x>SSL)  with  the 
inviscid  solution  governed  by 
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=  0 


with  pertaining  to  92v'/3y2  -  0  (i.e.,  vanishing  total  disturbance  shear) 

somewhere  within  the  Law  of  the  Wall  region12.  Once  the  v'(x,y)  field  is 
obtained,  the  attendant  streamwise  velocity  and  disturbance  shear  stress 
fields  may  be  found  from 


■  /X  if  d*  -  -  ^  if  O*  (»> 

-oo  _co  J 


An  important  and  useful  feature  of  this  approach  is  the  definition  of  an 
"effective  inviscid  wall"  position  (or  displacement  thickness)  that  emerges 
from  the  asymptotic  behavior  of  v'  far  from  the  wall8;  see  Figure  4.  As 
schematically  illustrated  in  Figure  4a,  this  is  defined  by  the  value  yweff 

where  the  "back  projection"  of  the  v.jnv  solution  vanishes,  this  projection 

being  the  (generally  non-linear)  solution  curve  obtained  by  inwardly- 
integrating  (18)  and  (12)  in  the  negative  y  direction  starting  at  5<^.  P  hy  si  - 

cal  1  y ,  yweff  thus  represents  the  total  mass  defect  height  due  to  the  shear 

stress  perturbation  field  and  hence  the  effective  wall  position  seen  by  the 
overlying  inviscid  middle  deck  disturbance  flow.  As  indicated  in  Figure  4b, 
this  concept  serves  to  couple  the  inner  and  middle  deck  solutions  in  a  direct 
physical ly-obvious  way  by  providing  the  non-singular  inner  equivalent  slip- 
flow  boundary  conditions  3p'/3y  (yeff)  =  vinv(yweff)  =  0  at  M^weff)  >  0  for 

the  middle  deck  solution  of  Eq.  (6).  In  practice,  this  approach  has  proven 
quite  useful  for  treating  a  variety  of  interaction  problems. 

C.  Approximate  Solution  by  Operational  Methods 

v 

An  analytical  solution  is  further  achieved  by  assuming  small  linearized 
disturbances  ahead  of,  behind  and  below  the  nonl i near  shock  jump  plus  an 
approximate  treatment30  of  the  detailed  shock  structure  within  the  boundary 
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SO.  Inger,  G.  R .,  "Shock  Wave  Penetration  and  Laterxil  Pressure  Gradient 

Effects  on  Transonic  Normal  Shock-Turbulent  Boundary  Layer  Interactions , 
AIAA  J.  15,  August  1977 ,  pp.  1179-1200 ,  plus  16 ,  May  1978 ,  p.  541 . 
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layer,  which  gives  reasonably  accurate  predictions  for  all  the  properties  of 
engineering  interest  when  >  1.05  (as  far  as  the  overall  interaction  proper¬ 
ties  are  concerned,  this  non-linear  shock  jump  provision  plus  the  various  non- 
uniform  viscous  flow  effects  within  the  boundary  layer  reduce  the  lower  Mach 
number  limit  otherwise  pertaining  to  the  linearized  supersonic  theory  in  pure¬ 
ly  inviscid  potential  uniform  flow).  As  described  in  detail  in  References  6 
and  10,  the  resulting  equations  can  be  solved  by  Fourier  transform  methods 
yielding  the  interactive  pressure  rise  and  displacement  thickness  growth 
inputs  to  the  above  extended  theory  of  the  inner-disturbance  sublayer.  The 
resulting  solution  contains  all  the  essential  physics  of  the  mixed  transonic 
viscous  interaction  field  for  non-separati ng  flows  including  the  upstream 
influence,  the  lateral  pressure  gradient  near  the  shock  and  the  onset  of 
incipient  separation;  numerous  detailed  comparisons  with  experiment31  have 
shown  that  it  gives  a  good  account  of  all  the  important  engineering  features 
of  the  interaction  over  a  wide  range  of  Mach-Reynolds  number  conditions. 

The  matching  of  the  outer  two  decks  with  the  inner  shear-disturbance  deck 
in  connection  with  the  Fourier  inversion  process  yields  the  determination  of 
the  upstream  influence  distance,  A  =  typical  numerical  results  for  it 

showing  the  important  typical  parametric  effects  of  Reynolds  number  and  shape 
factor  can  be  found  in  Ref.  6.  The  solution  further  yields  the  inner  deck 
displacement  thickness 
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and  the  interactive  skin  friction  relationship 
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where 
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31.  Inger ,  G.  R.,  "Application  of  a  Shock-Boundary  Layer  Interaction  Theory 
to  Transonic  Airfoil  Analysis ,"  AGARD  CP-291,  Colorado  Springs,  Sep  1980. 
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where  the  functions  H(T)  and  S(T)  are  given  in  Figure  5  and  represent  the  wall 
turbulence  effect  on  the  interactive  displacement  effect  and  skin  friction, 
respectively.  Figure  5,  in  fact,  is  a  central  result  of  the  gereral  non- 
asymptotic  theory,  providing  a  unified  account  of  the  entire  Reynolds  number 
range  in  terms  of  the  single  new  turbulent  interaction  parameter  T  (Eq.  23B) : 
it  ranges  from  the  limiting  behavior  of  negligible  wall  turbulence  effect  per¬ 
taining  to  Lighthill's  theory  at  T  +  0  (lower  Reynolds  numbers)  to  the  oppo¬ 
site  extreme  of  wall  turbulence-dominated  behavior  at  T»  1  pertaining  to  an 
asymptoti c- type  of  theory  at  vary  large  Reynolds  numbers  where  the  inner  deck 
thickness  and  its  disturbance  field  become  vanishingly  small.  The  relation¬ 
ship  between  these  two  heretofore-di sparate  theories  has  thus  been  explained 
and  established:  they  belong  at  opposite  extremes  of  a  general  non-asymptoti c 
theory  that  describes  the  transition  between  them.  Another  important  result 
emerging  from  Figure  5  is  that  the  asymptotic  trends  occurring  at  very  large 
Reynolds  numbers  cannot  be  extrapolated  down  to  ordinary  values;  doing  so  can 
yield  appreciable  error  in  the  inner  deck  properties  of  practical  interaction 
problems.  This  would  appear  to  explain  the  success  of  the  Lighthill  theory  in 
correlating  lower  Reynolds  number  turbulent  interactions  in  spite  of  the 
mathematical  rigor  of  asymptotic  theory:  the  former  is  simply  closer  to  the 
actual  physics  and  correctly  predicts  more  significant  interaction  effects  and 
scaling  under  the  decidedly  non-asymptoti c  conditions  involved.  By  the  same 
token,  the  extreme  approximation  involved  in  the  T  =  0  limit  significantly 
breaks  down  at  larger  Re^'s  pertaining  to  T»  1,  clearly  warranting  the  use 

of  the  present  theory  to  account  for  the  increasing  role  of  the  wall  turbu¬ 
lence  effect  on  the  interaction. 


An  important  and  useful  final  consequence  of  the  foregoing  analysis  is 
that  it  yields  an  explicit  analytical  criterion  for  the  onset  of  incipient 
separation  due  to  an  interactive  pressure  field;  setting  Tw  =  Tw  +  Tw  ~ 

Eq.  (22)  predicts  this  to  occur  when  o 
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where  it  is  re-emphasized  that  Cpw  here  is  the  local  interacti ve  pressure 

distribution.  Equation  (24)  bears  a  general  resemblance  to  a  Stratford- type9 
of  incipient  separation  relation  for  turbulent  flow,  except  that  the  present 
formula  contains  the  integrated  history  effect  along  the  interaction  whereas 
Stratford's  result  involves  purely  local  properties  of  Cp  and  dCp/dx.  It  is 

understood,  of  course,  that  the  present  theory  actually  breaks  down  approach¬ 
ing  such  separation  owing  to  the  linearization  assumptions  and  the  Van 
Driest/ Cebeci  wall  turbulence  model  used;  nevertheless,  Eq.  (24)  does  give  at 
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least  a  roughly-correct  indication  of  where  this  will  occur,  and  indeed  does 
so  without  containing  any  adjustable  empirical  constants. 

A  computer  program  has  been  constructed  to  carry  out  the  foregoing  solu¬ 
tion  method;  it  involves  the  middle-deck  disturbance  pressure  solution  coupled 
to  the  inner  deck  by  means  of  the  effective  wall  shift  (Eq.  21)  combined  with 
an  upstream  influence,  solution  subroutine.  The  corresponding  local  total 
interactive  displacement  thickness  growth  and  skin  friction  are  obtained  from 
Eqs.  (8)  and  (22),  respectively.  If  desired  the  attendant  boundary  layer 
shape  factor  change  along  the  interaction  may  also  be  calculated  as 
H  =  [S*Q  +  AS*(x)  J/0*(x)  with  6*  given  by  an  x-wise  integration  of  the  overall 

momentum  integral  equation  for  the  total  local  boundary  layer  since  p(x),  6* 
and  Cf  are  known.  The  incoming  turbulent  boundary  layer  is  treated  by  the 

compressible  version  of  a  universal  composite  Law  of  the  Wall  -  Law  of  the 
Wake  model  due  to  Walz32  that  not  only  has  a  convenient  analytical  form  (see 
Appendix)  but  also  provides  a  very  general  fundamental  description  of  this 
boundary  layer  in  terms  of  three  arbitrary  parameters:  preshock  Mach  number, 
boundary  layer  displacement  thickness  Reynolds  number,  and  the  incompressible 
shape  factor  H.  .  This  enables  us  to  account  for  the  important  influence  of 
nl 

the  upstream  flow  history  (pressure  gradient,  suction,  etc.)  on  the  inter¬ 
action. 

D.  Comparisons  of  the  Theory  with  Experiment 

Numerous  comparisons  of  the  present  theory  with  experimental  data  from 
both  wind  tunnel  and  free  flight  experiments  have  been  documented;  a  sample  is 
presented  in  Figure  6  to  illustrate  the  predicted  behavior  and  good  agreement 
for  the  interactive  pressure,  displacement  thickness  and  skin  friction  distri¬ 
butions  in  a  typical  non-separating  case.  Note  especially  in  Figure  6b  the 
lateral  pressure  gradient  that  occurs  in  the  vicinity  of  the  shock.  Regarding 
the  skin  friction  comparison  shown  in  Figure  6d,  we  note  that  the  "experimen¬ 
tal"  values  were  actually  inferred7  from  measured  velocity  profiles  along  the 
interaction  (hence  0*  and  6*)  by  means  of  the  2-D  momentum  integral  equation; 
considering  the  well-known  uncertainties  involved  in  their  experimental  set  up 
and  this  method,  combined  with  the  present  theory's  own  limitations,  the 
agreement  is  considered  good  as  regards  both  the  magnitude  and  shape  of  the  Cf 
curve. 

Figure  7  shows  a  comparison  of  the  predicted  local  shape  factor  change 
along  an  interaction  with  some  ONERA  measurements;  it  is  seen  that  the  theory 
nicely  captures  the  characteristic  local  peaking  of  H  near  the  shock  foot  as 
well  as  the  overall  behavior. 


32.  Walz ,  A .,  " Boundary  Layers  of  Flow  and  Temperature,"  M.I.T.  Press, 
Cambridge,  Mass.,  1969,  pp.  113. 
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A  particularly  interesting  comparison  with  some  very  recent  DFVLR-AVA  (GS) 
experiments33  is  shown  in  Figure  8,  where  we  compare  with  measurements  along 
the  non-separati ng  interaction  zones  on  two  supercritical  airfoils.  The 
experimental  Cf  values  were  inferred  from  streamwise  boundary  layer  profile 

surveys  by  means  of  the  Ludwi g-Ti 1 lman  relationship 


Cf  =  ,246(T*/Te) *796  Ree*"*268  e"1,561  Hi  (25) 


where  Hn-  is  given  by 

Hi  *  (H  -  .273Me2)/(l  .0  +  .1145Me2)  (26) 


p 

and  T*/Te  *  1  +  .14  M  for  y  =  1.40  and  P  »  .7  on  an  adiabatic  wall.  Theo- 

e  r 

retical  predictions  of  the  entire  airfoil  flow  field  were  also  made  with  a 
composite  inviscid  transonic  -  turbulent  boundary  layer  -  shock  interaction 
numerical  scheme  in  which  the  present  theory  was  used  as  a  local  interactive 
module  astride  the  inviscid  shock  location  (see  Reference  33  for  details).  It 
is  seen  that  the  theory  yields  an  excellent  prediction  of  the  local  skin  fric¬ 
tion  values  upstream  and  slightly  downstream  of  the  shock  (including  the  mini¬ 
mum  value)  and  a  good  qualitative  account  of  the  overall  shape  of  the  stream- 
wise  distribution.  Weil  downstream  of  the  shock,  the  theory  slightly  over¬ 
estimates  the  post-shock  Cf  recovery  inferred  by  the  Ludwi  g-Ti 1 lman  relation¬ 
ship.  Nevertheless,  in  view  of  the  combined  limitations  of  the  experimental 
Ludwi  g-Ti 1 lman  method  and  the  present  small  disturbance  theory,  the  overall 
agreement  with  the  data  is  regarded  as  good. 


III.  EXTENSION  OF  THE  LOCAL  INTERACTION  THEORY 
TO  SPINNING  AXISYMMETRIC  BODIES 

A.  Fundamental  Assumptions 

For  the  high  Reynolds  number  non-separating  flows  at  zero  angle  of  attack 
which  we  are  considering  here,  the  extension  of  the  foregoing  interaction 
theory  to  spinning  bodies  of  revolution  is  based  on  the  following  general 
physical  assumptions.  (1)  The  three  dimensional  relief  and  spin  effects, 
while  altering  certain  details  and  quantitative  values  of  the  parameters 


33.  Nandanan ,  M. ,  St anew sky,  E.  and  Ingev,  G.  R.,  "A  Computational  Procedure 
for  Transonic  Airfoil  Flow  Including  a  Special  Solution  for  Shock- 
Boundary  Layer  Interaction,  "  AIAA  Journal  Vol.  19,  dec.  81,  pp.  1540-46. 
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involved,  do  not  fundamentally  change  the  triple  deck  interaction  structure 
illustrated  in  Figure  3.  Some  detailed  analysis  has  shown  that  this  might  not 
be  true  only  for  very  slender  "needle-like"  bodies  with  unrealistically  -  high 
spin  rates.  (2)  Being  a  small  scale  local  event  dominated  by  rapid  changes 
in  a  direction  normal  to  the  inviscid  shock,  the  physics  of  the  interaction  is 
still  essentially  "two-dimensional"  even  on  a  spinning  axi symmetric  body;  the 
significant  effects  of  the  three  dimensional  geometry  and  spin  enter  primarily 
from  their  influence  on  the  surrounding  large  scale  inviscid  flow,  shock  loca¬ 
tion  and  incoming  boundary  layer  properties  ahead  of  the  interaction.  (3) 
Excluding  exceedingly  slender  bodies  and/or  very  low  Reynolds  number  condi¬ 
tions,  both  the  local  turbulent  boundary  layer  thickness  6Q  and  streamwise 

scale  of  the  interaction  zone  (5-7  SQ)  are  small  enough  compared  with  the 

local  transverse  and  longitudinal  body  radii  of  curvature  to  permit  neglect  of 
the  explicit  longitudinal  and  viscous  transverse  curvature  terms  in  the 
governing  flow  equations  and  to  permit  retention  of  the  parallel  shear  flow 
approximation  for  the  incoming  boundary  layer  as  far  as  the  interaction  solu¬ 
tion  is  concerned.  (4)  The  body  surface  and  the  overlying  boundary  layer 
flow  are  adiabatic. 

We  now  proceed  to  examine  the  basic  relations  upon  which  the  interaction 
solution  rests,  to  bring  out  specifically  the  influence  of  axi symmetric  flow 
geometry  and  spin  effects. 

B.  The  Inviscid-Disturbance  (Middle  and  Outer  Deck)  Regions 

1.  Some  Fundamental  Relationships.  For  the  case  of  bodies  at  zero  angle 
of  attach  the  general  equations  governing  the  attendant  circumferentially  - 
invariant  (3/3<j>  =  0)  flow  field  are  as  follows  when  expressed  in  the  body- 
oriented  coordinates  shown  in  Figure  9  (such  cordi nates  have  been  chosen  here 
because  of  their  convenience  in  treating  the  viscous-interaction  aspects  of 
the  problem) : 
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CIRCUMFERENTIAL  MOMENTUM 
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ENERGY  CONSERVATION 


(31) 


where  a2  =  9p/9p  =  yp/p,  r  =  rB(x)  +  ycos  Og  and  where  for  the  moment  we  have 

included  the  viscous  terms  to  the  boundary  layer-order  of  approximation. 
These  equations  are  of  course  supplemented  by  the  thermal  equation  of  state. 
Note  also  that  Eq.  (30)  governs  the  spin-induced  circumferential  (cross)  flow 
caused  by  viscous  shear  in  the  boundary  layer. 

We  now  proceed  to  obtain  an  important  gas  dynamic  relationship  involving 
the  velocity  field  components  as  follows.  First  we  note  from  continuity  (27) 
that 
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which  when  substituted  into  the  energy  Equation  (31)  yields  the  convective 
derivative  of  pressure 


(33) 


Now  since  the  normal  and  circumferential-flow  Mach  numbers  v/a  and  w/a, 
respectively,  are  very  small  in  the  present  problem  we  can  hereafter  neglect 
their  squares  compared  with  unity  with  excellent  approximation;  doing  this 
while  substituting  9p/9x  from  Equation  (28)  plus  9p/9y  from  Equation  (29)  into 
the  LHS  of  Equation  (33),  we  obtain  after  some  algebra  the  following  equation: 
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where  5  5  9u/3y  -  9v/9x  is  the  vorticity.  Equation  (34)  is  the  generalized 
version  of  the  classical  flow  equation  of  compressible  inviscid  flow  which 
here  explicitly  accounts  for  axisymmetric  flow  effects  and  the  possible  pres¬ 
ence  of  vorticity;  it  thus  pertains  to  flow  inside  as  well  as  outside  the 
boundary  layer. 

A  second  companion  equation,  of  particular  use  in  interaction  problems,  is 
the  generalized  streamline  slope  relationship  obtained  by  substituting  3u/3x 
from  Equation  (28)  into  the  RHS  of  Equation  (33)  and  rearranging  the  result 
into  an  equation  governing  v/u;  again  dropping  terms  of  order  (v/a)2  and 
(w/a)2,  we  get 
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Equation  (35)  enables  one  to  calculate  v/u  across  the  flow  by  appropriate 
integration  of  the  pressure  and  viscous  stress-gradient  effects  including  the 
3-D  relief  effect  from  the  last  LHS  term. 

A  third  important  relationship  governs  the  pressure  field  in  a  vortical 
but  otherwise  negligible  shear  stress  region  of  flow.  This  we  obtain  by  drop¬ 
ping  the  explicit  stress  terms  in  solving  Equation  (35)  for  3p/3x  and  then 
differentiating  the  result  to  obtain,  after  re-use  of  Eq.  (35),  an  expression 
for  32p/3x2;  this  expression  is  then  subtracted  from  the  result  of  differenti¬ 
ating  Eq.  (29)  w.r.t.y.  After  much  algebra  involving  the  repeated  use  of  Eqs. 
(28),  (29)  and  (35)  plus  the  neglect  of  all  v2/a2  terms  we  ultimately  obtain 
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This  equation  constitutes  a  very  generalized  form  of  Lighthill's  famous  pres¬ 
sure  field  equation  extended  here  to  include  arbitrary  non-parallel 
(v,3u/3x  *  0)  vortical  shear  flows  on  2-D  or  spinning*  axisymmetric  bodies. 

Equations  (34),  (35)  and  (36)  comprise  a  trio  of  basic  relations  which 
govern  the  total  interactive  flow  field  in  the  combined  outer-middle  deck 
regions.  We  now  specialize  them  to  apply  to  the  interactive  disturbance  field 
by  imagining  that  the  flow  consists  of  an  incoming  turbulent  boundary  layer 
plus  an  overlying  inviscid  transonic  flow  (subscript  "o")  which  is  slightly 
perturbed  by  the  interaction  caused  from  a  weak  nearly-normal  shock  standing 
in  the  outer  flow;  we  thus  write  u  =  UQ  +  u'(x,y) ,  v  =  vQ  +  v',  w  =  wQ  +  w-* , 

p  =  PQ  +  p'  and  p  =  pQ  +  p'  where  the  primed  quantities  denote  these  inter¬ 
active  perturbations.  Substituting  these  perturbation  expressions  into  Eq. 
(34),  neglecting  all  shear  stress  disturbances  by  virtue  of  the  "frozen  turbu¬ 
lence"  approximation,  assuming  there  is  negligible  correlation  between  turbu¬ 
lent  fluctuations  and  these  perturbations,  and  subtracting  out  the  undisturbed 
flow  equations  and  linearizing  the  resulting  small  perturbation-fl ow  equations 
except  for  the  appropriate  transonic  terms,  we  ultimately  get  the  following 
disturbance  flow  velocity  equation: 
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where  for  zero  angle  of  attack  w'  =  0  since  the  shock  position  is  independent 
of  4s  and  hence  the  resulting  interaction  field  cannot  induce  any  circumferen¬ 
tial  cross  flow  regardless  of  spin.  Now  in  a  typical  transonic  flow  field 
outside  the  inner  deck,  the  terms  of  the  last  square  bracket  on  the  LHS  of 
(37A)  are  very  small  compared  with  the  remaining  ones  and  hence  may  be  neglec¬ 
ted;  likewise,  the  terms  proportional  to  v0/a0  are  negligiably  small  compared 

with  v'/aQ.  Consequently  to  the  leading  order  of  approximation,  Eq.  37A  may 

be  simplified  to 


*The  spin  here  is  implicit,  of  course,  since  all  the  explicit  spin  terms  are 
of  negligible  order  w^/a^  and  hence  have  been  dropped. 
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which  is  seen  to  be  a  generalized  transonic  small  disturbance  equation  for 
either  a  2-D  or  axi- symmetric  rotational  perturbation  field  in  a  parallel 
background  flow  M0(y).  Proceeding  in  a  similar  manner,  the  application  of  Eq. 

(35)  to  the  small  disturbance  field  yields  the  following  approximate  linear¬ 
ized  expression  governing  the  perturbation  streamline  slope  vVUQ: 
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while  the  comparable  perturbation  version  of  (36)  is  found  to  be 
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where  the  RHS  here  is  evidently  a  kind  of  accoustical  source-disturbance 
associated  uniquely  with  the  3-D  geometry  in  a  compressible  flow. 

Equations  (37B),  (38)  and  (39)  are  the  three  basic  equations  governing  the 
transonic  interactive- small  disturbance  field  (u',  v^,  p^)  in  both  the  outer 
and  middle  decks  on  either  2-D  or  spinning  axi-symmetric  bodies  at  zero  angle 
of  attack.  They  offer  the  advantage  that  the  three-dimensional  effects  appear 
explicitly  and  hence  their  influence  can  be  assessed  directly. 

2.  Influence  of  3-D  Effects.  In  preparation  for  examing  these  interac¬ 
tion  equations,  we  first  evaluate  certain  terms  that  appear  in  the  coeffi¬ 
cients.  Using  r  =  rg(x)  +  y  cos0g  with  rg  the  body  radius  and 

tan0g  =  drg/dx,  we  obtain 
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Now  consistent  with  our  previous  neglect  of  longitudinal  curvature  terms  in 
the  governing  equations,  we  may  also  neglect  such  terms  here  (i.e.,  those 
~  d0g/dx);  thus  r"1  8r/3x  has  at  most  the  "near  body"  value  (40B)  while  being 

considered  negligible  in  the  far  field  according  to  (40C).  Using  Eqs.  (40) 
and  (41)  we  further  obtain  the  estimates 

|_(lr/3X)  _  (irpxjz  «  _2r-2  sin2  0B  near  body  (42A) 

=  0  in  far  field  (42B) 

and 
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where  we  note  that  for  body  regions  where  0g«  1  the  values  given  by  Eqs.  (42) 
and  (43)  are  very  small  indeed. 

Now  consider  the  transonic  small  perturbation  Equation  (37B),  which  we 

note  includes  shock  jump  effects  via  the  non-linear  term  (y  +  1)  MQ2  y”  "ax-* 
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Within  the  near  field  region  y  <  6  pertaining  to  the  middle  deck,  we  have  the 
order  of  magnitude  estimates  that 
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which  in  turn  show  that  for  slender  bodies  the  explicit  3-D  term  in  Eq.  37A  is 
of  order  6Q/r g  smaller  than  the  remaining  terms  and  hence  can  be  neglected  in 

the  leading  approximation  for  high  Reynolds  number  turbulent  boundary  layer 
flows  where  typically  SQ/rB  <  .05  -  .10.  This  implies  that  the  important 

middle  deck  region  is  dominated  by  the  two-dimensional  aspects  of  the  short- 
scale  interaction.  Turning  to  the  intermediate  field  region  y  ~  rD  ~  r  within 

D 

the  inner  portion  of  the  outer  deck,  we  have  the  order  of  magnitude  estimates 
that 
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where  Xinter  is  the  streamwise  scale  of  the  interaction.  Since 
5  <  Xinter/6Q  <  10  is  much  smaller  than  rB/6Q  in  practical  applications,  we 

see  that  the  term  (45B)  is  much  larger  than  the  other  two,  again  implying  as 
far  as  Eq.  37 B  is  concerned  that  the  local  physics  is  predominantly  two- 
dimensional  when  we  include  the  local  transonic  shock  jumps  (as  indeed  we  do) 
in  the  treatment  of  the  outer  deck.  In  the  far  field  region  y  ~  r  >  Xinter  in 
the  outer  portion  of  the  outer  deck,  we  estimate  that 
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which  for  r  >  Xinter  are  all  of  the  same  magnitude  (except  for  u'9r/9x) . 
Consequently,  in  this  region  the  explicit  radial  spreading  term  v'9r/9y  in  Eq. 
37B  now  becomes  of  the  same  order  as  the  remaining  2-D  terms  in  governing  the 
inviscid  disturbance  field;  as  is  well-known  the  axisymmetric  geometry  does 
indeed  influence  the  asymptotic  behavior  in  the  far  field  by  increasing  the 
disturbance  decay  rate  over  the  2-D  result.  In  the  present  application,  how¬ 
ever,  where  the  finer  details  of  this  weak  far  field  behavior  are  not  of 
interest  nor  of  too  much  consequence  as  regards  the  aerodynamic  forces  and 
moments,  we  shall  still  neglect  the  explicit  3-D  term  in  the  outer  deck  solu¬ 
tion  of  Eq.  37B.  Indeed,  any  modest  errors  involved  in  this  approximation  are 
counteracted  by  the  fact  that  the  resulting  " 1 ocal ly-2-D"  interaction  is 
centered  about  the  correct  shock  position  from  the  axisymmetric  body  inviscid 
sol uti on. 

A  similar  examination  of  the  disturbance  streamline  relationship  (38) 
using  the  order  of  magnitude  estimates  (40)  -  (46)  yields  the  same  conclu¬ 
sions:  the  explicit  3-D  geometry  terms  are  an  order  of  magnitude  smaller  than 
the  remaining  interaction-effect  terms  in  both  the  middle  deck  and  near-field 
portion  of  the  outer  deck,  and  only  in  the  far- field  region  y  ~  r  >  Xinter  of 
the  outer  deck  (where  in  fact  v'  is  quite  weak)  is  the  3-D  spreading  term  in 
(38)  important.  An  examination  of  the  pressure  perturbation  Eq.  (39)  again 
using  (40)  -  (46)  shows  that  the  entire  RHS  is  negligibly  small  for  slender 
bodies  regardless  of  the  y  scale,  the  only  explicit  3-D  term  of  possible 
significance  being  the  radial  spreading  effect  in  the  coefficient  of  9p/9y  on 
the  LHS.  By  the  aforementioned  arguments,  this  latter  effect  is  important 
only  in  the  far-field  weak-di sturbance  region  of  the  outer  deck. 

To  summarize,  it  has  been  shown  that  for  high  Reynolds  number  transonic 
flows  on  slender  axisymmetric  bodies  of  practical  interest  the  dominant  shock¬ 
boundary  layer  interaction  physics  is  locally  two-dimensional  except  in  the 
weak  far-field  region  of  the  outer  deck  where  the  explicit  3-D  spreading- 
effect  terms  are  influential  in  determining  the  asymptotic  disturbance  decay. 
Consequentially,  to  a  good  engineering  approximation  in  treating  the  local 
shock-boundary  layer  interaction  effects  on  projectile  flow  fields  of  prac¬ 
tical  interest,  the  2-D  form  of  the  interaction  equations  (as  solved  already 

in  Ref.  6)  may  be  used  on  slender  bodies  of  revolution  as  well.  Indeed, 

unless  the  finer  details  of  the  far  field  are  of  interest  (as  for  example  in 

sonic  boom  determination)  the  resulting  minor  errors  are  alleviated  by  the 

fact  that  both  the  coeffi cients  in  these  2-D  equations  and  the  far-field 
inviscid  shock  location  fully  include  the  axi-symmetric  geometry  effects.  In 
a  similar  vein,  it  was  shown  that  in  practice  body  spin  adds  no  explicit  new 
terms  to  the  governing  equations  for  flight  at  zero  angle  of  attack;  the  only 
effects  (which  we  fully  account  for)  are  those  implicit  through  the  spin 
effect  on  the  boundary  layer  properties  that  appear  as  coefficients  in  the 
disturbance  equations. 
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C.  The  Inner  Shear-Disturbance  Deck 


The  extreme  thinness  of  this  sublayer  immediately  justifies  the  neglect 
therein  of  transverse  curvature  effects  [i.e.,  r  «  rg(x)],  lateral  pressure 

gradients  and  temperature  variations  above  the  adiabatic  surface.  We  may  also 
use  the  boundary  layer-approximation  to  the  shear  stress  field.  Taking  the 
layer  to  lie  within  the  Law  of  the  Wall  portion  of  the  incoming  turbulent 
boundary  layer  as  in  the  2-D  case,  then,  the  appropriate  governing  equations 
of  the  total  interactive  flow  on  a  spinning  axi- symmetric  body  at  zero  angle 
of  attack  become: 
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where  Eq.  (49)  governs  the  distribution  of  the  spin-induced  cross  flow  w0(y) 
in  the  incoming  non-i nteracted  boundary  layer.  Here  p  and  v  are  the  den- 
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sity  and  kinematic  laminar  viscosity  coefficient,  respectively,  based  on  adia¬ 
batic  wall  conditions,  £jxx  and  ape  the  streamwise  and  circumferential 

eddy  viscosity  functions,  respectively,  and  the  known  pressure  pw(x)  is  given 

by  the  solution  along  the  bottom  of  the  middle  deck.  In  the  light  of  the 
results  from  previous  boundary  layer  studies  on  spinning  projectiles1,  it  will 
be  assumed  here  (excluding  abnormal ly-hi gh  spin  rates  and/or  extremely  slender 
bodies)  that  neither  the  axi-symmetri c  geometry  nor  spin  fundamentally  alters 
the  turbulent  structure  near  the  wall;  accordingly  it  is  sufficient  for  the 
present  purposes  to  adopt  the  isotropic  VanDriest-Cebeci  eddy  viscosity  model 


eTxx~  eT 4><|> 


-y/LD  2 

k2y2  (1  -  e  ) 


+ 


aw 

(— )2 
'  3y  1 


(51A) 


with 


32 


( 5 1 B) 


Ln  =  A  v  (T 

H  n  '  u 


w 


’TOT 


/pwo) 


-1/2 


w. 


=  /t 


TOT 


,  +  T  2 

“x  '"I 


(51C) 


where  tw  =  uwq  (3u/3y)w,  tw  =  uwQ  (3w/3y)w,  k  =  .41  and  A  =  26. 

X  <() 

Consider  first  the  undisturbed  zero  pressure  gradient  flow  in  the  incoming 
turbulent  boundary  layer  close  to  the  wall.  Although  the  presence  of  cross 
flow  wQ(y)  due  to  spin  alters  the  usual  constancy  of  total  shear  across  the 

Law  of  the  Wall  region  (as  can  be  seen  from  Eqs.  (48)  and  (49)  with 

wo  ~  wo(°)  +  V<°>y.  wo(0)  -  n  Tg  in  the  left  hand  sides),  further  analysis 

shows  this  effect  to  be  negligibly  small  for  the  high  Reynolds  number-spinning 
slender  body  flows  of  practical  interest.  Hence  to  a  good  approximation  the 
undisturbed  profiles  U0(y) ,  w0(y)  obey  the  following  relations  obtained  by 

integrating  Eqs.  48  and  49  after  neglecting  their  left  hand  sides: 
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where  eT  and  eT  are  given  by  Eqs.  (51)  with  U(y)  =  U0(y)  and  where  the 
xxo  *+0 

wall  stress  ratio  a  =  (t  /t  )  is  presumed  a  known  function  of  wg  =  wo(0)  = 
.  <t>  x 

n  Tg  which  vanishes  as  w^  o. 


We  now  imagine  the  shock-boundary  layer  interaction  process  to  induce 
small  perturbations  u'  =  u  -  UQ,  v'  =  v,  p'  =  p  -  pQ  without,  however,  affect¬ 
ing  the  cross  flow  (w'  *  o,  see  above).  Then  restricting  attention  to  the 
non-separating  case  by  linearizing  these  disturbances  in  the  first  approxima¬ 
tion  and  utilizing  the  parallel  shear  flow  approximation  (Vo  =  8Uo/8x  -  0)  for 
the  undisturbed  flow  consistent  with  Eqs.  52  and  53,  Eqs.  47  and  48  yield  the 
following  disturbance  flow  relations: 
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is  a  measure  of  the  local  spin-induced  cross  flow  angle,  which  evidently 
influences  the  streamwise  turbulent  stress  perturbation  field  near  the  sur¬ 
face.  Since  a  for  practical  spinning  projectiles  is  typically  quite  small 
(a  <  .05) 1,  however,  it  is  seen  that  its  effect  (~  a2)  on  the  inner  disturb¬ 
ance  solution  is  negligible:  Eqs.  (54)  thru  (56)  thus  reduce  to  their  non¬ 
spinning  form. 

We  now  assess  the  explicit  influence  of  axi symmetry  on  the  inner  solution, 
which  as  can  be  seen  enters  only  via  the  disturbance  continuity  equation  (54). 
That  this  influence  is  negligible  for  slender  bodies  over  the  relatively  short 
streamwise  scale  Xinter  and  very  small  inner  deck  height  seal  e  5sl  can  be 
established  by  comparing  the  following  order  of  magnitude  estimates  of  the 
various  terms  in  Eq.  (54): 
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Clearly,  when  Xinter/rB  and  <So/rB  are  small,  the  last  terms  in  Eq.  (58)  are 

negligible  for  slender  bodies  and  Eqs.  (54)  thru  (56)  thus  reduce  to  their 
two-dimensional  form  (with  however  the  coefficients  based  on  spinning  axisym- 
metric  solution  values);  like  the  outer  and  middle  decks,  the  dominant  inter¬ 
action  physics  of  the  inner  deck  is  locally  two-dimensional. 
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D.  The  Overall  Displacement  Thickness 


The  integrated  continuity  equation  across  the  boundary  layer  region  on  an 
axi symmetric  body  yields  the  following  general  expression  for  the  total  dis¬ 
placement  thickness  6*  =  /®(1  -  _£g_)  dy  assuming  isentropic  inviscid  flow 
along  the  boundary  layer  edge:  e 
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where  ve/ue  is  the  streamline  slope  at  the  edge  and  where  the  transverse 
curvature  effect  in  the  integrand  has  been  neglected  for  S/rg  «  1.  The 

corresponding  interactive  small  perturbation  part  of  Eq.  (59),  using  the  2-D 
version  of  Eq.  (38)  integrated  across  the  boundary  layer,  gives  the  local 
interactive  displacement  thickness  increase  as 
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Now  the  first  bracket  on  the  right  side  of  Eq.  (60B)  contains  two  contribu¬ 
tions  from  the  interactive  pressure  field;  the  first  involves  a  Mach  number 
profile  integral  across  the  middle  deck  plus  the  inner  deck  contribution  from 
the  non-zero  layer  limit  yw  eff  as  previously  explained,  whereas  the  second  is 

a  purely  inviscid  mass  flow  effect.  The  last  term  on  the  right  constitutes  a 
small  explicit  contribution  from  the  axi symmetric  spreading  effect  which  is 
negligible  in  our  applications  because  it  involves  the  product  of  two  terms 
which  are  both  very  small  in  the  interaction  regions  on  a  slender  body  at  high 
Reynolds  numbers;  discarding  it  thus  leaves  the  two-dimensional  form  of  Eq. 
(60B)  in  which  axisymmetry  and  spin  effects  are  implicitly  accounted  for  in 
the  coefficient  values  of  6o,  6o*,  Mei,  p'  and  Mo(y) . 
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IV.  COMPOSITE  VISCOUS-INVISCID  FLOW  FIELD  MODEL 


A  composite  model  of  the  complete  flow  field  around  spinning  projectiles 
at  supercritical  flight  speeds  has  been  constructed  which  utilizes  the  forego¬ 
ing  shock-boundary  layer  interaction  theory  as  a  locally-imbedded  interactive 
module  astride  each  body  shock  location  (see  Figures  1,  3).  We  now  describe 
the  other  components  of  this  model  and  how  they  are  coupled  with  the  interac¬ 
tion  zone  solution. 

A.  Inviscid  Solution 


The  inviscid  flow  calculations  are  treated  by  the  method  developed  by 
Reklis,  Sturek ,  and  Bailey5  involving  the  numerical  solution  of  the  transonic 
small  disturbance  equation  for  the  velocity  perturbation  potential  $  which  in 
cylindrical  coordinates  is 


[(i-M^)  -  m£  (y+1)$x]$xx  +  *rr 
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This  equation  is  a  second  order  nonlinear  partial  differential  equation  of 
mixed  el  1 i pti c/hyperbol i c  type;  the  type  changes  to  match  the  physical  differ¬ 
ences  between  regions  of  subsonic  and  supersonic  flow. 

Equation  (61)  generally  yields  adequate  predictions  of  the  inviscid  flow 
about  a  projectile  shape  such  as  that  studied  here.  Certain  regions  of  the 
flow  require  some  "modeling",  however.  The  wake  is  simulated  by  a  necked-down 
extension  from  the  boattail  base  smoothly- fared  into  a  cylindrical  sting  over 
a  distance  of  two  calibers;  the  base  flow  is  thus  modeled  as  an  extended 
sting.  A  review  of  free- flight  shadowgraphs  for  projectile  shapes  at  transon¬ 
ic  speeds  does  show  the  wake  flow  to  follow  near  the  boattail  angle  for  a 
distance  of  one  to  three  calibers  before  turning  parallel  to  the  flow  direc¬ 
tion,  but  of  course  the  wake  vorticity  which  physically  exists  is  not  account¬ 
ed  for.  Moreover,  this  obviously  ignores  the  wake  momentum  defect  and  hence 
the  base  contribution  to  the  overall  drag,  as  well  as  slightly  over-estimating 
the  pressure  rise  and  hence  the  boundary  layer  thickness  growth  and  skin  fric¬ 
tion  drop  near  the  end  of  the  boattail  owing  to  neglect  of  the  base  pressure 
upstream  influence.  However,  for  the  present  purposes  where  drag  prediction 
is  not  of  primary  interest,  these  shortcomings  are  not  deemed  too  significant. 

In  order  to  develop  a  "conservative"  algorithm  to  solve  this  equation 
special  care  must  be  taken  at  transitions  between  subsonic  and  supersonic 
flow.  Non-conservative  forms  of  the  algorithm,  however,  often  give  better 
agreement  with  experiment  because  the  breakdown  in  conservation  approximates 
the  effect  of  the  shock  obliquity  resulting  from  shock-boundary  layer  interac¬ 
tion  (see  below).  The  algorithm  used  here  is  therefore  a  non-conservati ve  one 
for  purposes  of  a  first  order  composite  flow  model.  Consistently,  first  order 
boundary  condition  relations  have  been  used  with  120  streamwi se  grid  points 
along  the  body  length  with  higher  concentration  in  the  corner/shock  inter¬ 
action  regions. 
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B.  Boundary  Layer  Theory 


Turbulent  boundary  layer  flow  computations  were  made  by  an  extension  of 
the  method  originated  by  Dwyer  and  Sanders1’34.  In  this  technique  the  bound¬ 
ary  layer  equations  for  the  conservation  of  mass,  momentum,  and  energy  are 

solved  with  an  implicit  finite  difference  technique.  The  solution  begins  with 
the  development  of  an  approximate  boundary  layer  profile  at  the  tip.  The 

solution  is  then  marched  along  the  body  from  nose  to  tail.  At  each  step  along 
the  way  a  two  point  boundary  value  system  is  solved  with  conditions  given  at 
the  body  surface  and  at  the  boundary  layer  edge.  The  possibility  of  body  spin 
is  accounted  for  and  care  is  taken  in  setting  up  the  difference  equations 
which  are  solved  to  maintain  differencing  in  the  correct  direction  to  maintain 
stability.  The  solution  is  started  at  the  body  nose  tip  with  the  assumption 
that  it  is  a  cone  and  that  the  boundary  layer  solution  is  self  similar  along 
rays  through  the  cone  tip.  This  assumption  is  strictly  true  only  if  the  flow 
at  the  tip  is  supersonic.  Errors  will  damp  out  as  the  solution  proceeds  down 
the  body,  however,  and  even  when  the  tip  is  emersed  in  subsonic  flow  the 

starting  conditions  obtained  from  this  method  have  been  satisfactory.  Turbu¬ 
lence  is  accounted  for  by  use  of  a  turbulent  shear  stress  model  with  Van 
Driest  damping.  This  model  has  proved  suitable  for  use  in  cases  of  supersonic 
flow  and  is  carried  over  directly  to  the  transonic  regime.  Although  this 
method  of  solving  the  boundary  layer  equations  has  been  used  for  cases  of 
supersonic  flow  over  projectiles1  in  the  past,  this  was  the  first  use  of 
Dwyer's  method  for  problems  in  transonic  flow. 

C.  Imbedded  Local  Interaction  Solutions 


The  presence  of  a  local  compression  shock  abruptly  terminating  the  local 
pocket  of  supersonic  flow  following  both  body  corners  is  identified  and  loca¬ 
ted  by  the  inviscid  solution  code  and  used  to  call  the  interaction  solution  as 
an  interuptive  subroutine  centered  about  each  such  shock  location,  using  the 
appropriate  local  pre- shock  Mach  number  from  the  code  as  a  driving  input  to 
the  solution.  When  inserting  the  interaction  we  account  for  the  fact  that  the 
non-conservative  inviscid  code  gives  shock  strengths  slightly  less  than  the 
full  Ranki ne-Hugoni ot  normal  shock  jump  assumed  in  the  interaction  solution  by 
correcting  the  pre-shock  Mach  number  M}  used  in  the  latter  solution  to  the 

value  pertaining  to  a  slightly  oblique  attached  shock  that  produces  maximum 
post- shock  flow  deflection;  this  method,  which  is  strongly  supported  by  a 
correlation  of  experimental  evidence  by  many  investigators  (see  e.g..  Ref.  33 
and  Fi gure  11) ,  yields 


Mi  eff  “  Mi  si n(90  -  37.8  /Mx  -  1  ) 


(62) 


34.  Dwyer ,  H.  A .,  and  Sanders ,  E.  7?.,  "Magnus  Forces  on  Spinning  Supersonic 
Cones.  Part  I.  The  Boundary  Layer  >  AIAA  Journal  14,  April  1976, 
p.  498. 
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The  additional  required  inputs  of  the  incoming  boundary  layer  displacement 
thickness  and  shape  factor  are  given  by  the  aforementioned  turbulent  boundary 
layer  code  that  has  been  simultaneously  calculating  along  the  body.  The 
interaction  subroutine  replaces  the  boundary  layer  code  over  the  range  of  the 
interaction  with  a  complete  description  (if  desired)  of  the  local  small  scale 
interaction  wall  pressure,  displacement  thickness  and  skin  friction  distribu¬ 
tion  (such  as  those  shown  in  Figure  6),  plus  the  final  post-interaction  values 
of  6*,  Cf  and  the  (subsonic)  inviscid  edge  conditions  needed  to  restart  the 
turbulent  boundary  layer  code  downstream.  This  method  of  introducing  the 
interactive  solution  allows  us  to  account  not  only  for  the  rapid  displacement 
thickness  growth  in  each  interaction  but  also  for  the  attendant  interactive 
distortion  of  both  the  skin  friction  and  profile  shape.  Moreover,  the  import¬ 
ant  influence  of  these  changes  on  the  subsequent  boundary  layer  development 
downstream  is  included  in  appropriate  post-interaction  reinitialization  of  the 
turbulent  boundary  layer  calculation.  Consistent  with  the  turbulent  boundary 
layer  model  employed  within  the  interaction  theory  to  fit  the  pre-interaction 
flow  conditions,  this  reinitialization  is  carried  out  at  the  first  post- 
interactive  streamwise  station  by  means  of  the  compressible  version  of  Walz's 
general  composite  Law  of  the  Wall  -  Law  of  the  Wake  velocity  profile  model 
(Appendix  A)  with  the  profile  parameters  chosen  to  match  the  post-interactive 
values  of  6*,  Cf  and  Me  given  by  the  local  interaction  solution,  leaving  the 

wake  function  and  shape  factor  to  take  their  consequent  highly  non-equilibrium 
post-shock  values.  Given  this  profile  fit,  the  turbulent  boundary  layer  code 
then  marches  downstream. 

It  should  be  noted  that  since  the  interaction  zones  on  practical  artillery 
shell  result  from  the  rapid  overexpansi on-recompressi on  regions  following  the 
ogive-cylinder  and  ogi ve-cyl inder-boattai 1  junctions  (see  Figures  1  and  10) 
the  aforementioned  general  turbulent  boundary  layer  profile  model  in  the 
interaction  module  is  required  to  account  for  the  important  effect  of  the 
rapid  shape  factor  variation  along  these  regions.  This  implies  that  a  suffi¬ 
ciently  large  number  of  small  boundary  layer  x-steps  be  concentrated  in  such 
regions.  On  the  other  hand,  these  corner-zones  simplify  the  flow  problem  as 
well:  they  fix  the  shock  location  nearly  independently  of  the  viscous  dis¬ 

placement-thickness  effects.  An  illustration  of  this  is  shown  in  Figure  12, 
where  typical  Navier-Stokes  code  pressure  distribution  predictions35  at  dif¬ 
ferent  Reynolds  numbers  are  shown  to  be  close  to  the  inviscid  result  as  corro¬ 
borated  by  experiment.  Thus,  in  contrast  to  the  case  of  a  supercri ti cal  wing 
where  the  shock  location  and  strength  are  unknown  a  priori  and  significantly 
influenced  by  boundary  layer  displacement  effects  (thereby  requiring  a  global 
vi scous-inviscid  interaction  calculation36),  the  present  problem  to  a  very 


35.  Nietubicz,  C.  J .,  Pulliam,  T.  H.,  and  Steger,  J.  L.,  "Numerical  Solution 
of  the  Azimuthal-Invariant  Thin-Layer  Navier-Stokes  Equations  ARBRL-TR- 
02227,  March  1980.  AIAA  Paper  79-0010,  January  1979. 

36.  Stanewsky,  E.,  Nanandan,  N.,  and  Inger,  G.  R.  ffroc.  AGARD  Symposium  on 
"Computation  of  Viscous-Inviscid  Interactions >  AGARD  CP-291,  Colorado 
Spring,  September  1980. 
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good  approximation  can  take  the  shock  location  given  by  the  purely  inviscid 
solutions  without  any  global  iteration. 


V.  DISCUSSION  OF  THEORETICAL  RESULTS  AND  COMPARISONS  WITH  EXPERIMENT 

The  present  composite  vi scous-i nvi scid  interaction  theory  has  been  devel¬ 
oped  into  a  global  flow  field  prediction  computer  code  that  is  fully  opera¬ 
tional  at  BRL.  It  has  been  exercised  on  a  number  of  supercritical  projectile 
cases  at  zero  angle  of  attack,  for  which  experimental  data  has  also  been 
obtained;  following  a  brief  description  of  these  experiments,  we  will  present 
comparisons  below.  We  will  also  compare  with  predictions  based  on  a  thin- 

layer  Navier-Stokes  numerical  code  used  by  BRL.  Finally,  we  will  show  some 
interaction-effect  predictions  in  the  presence  of  projectile  spin. 

A.  Experimental  Data. 

To  provide  a  sound  validation  base  for  the  development  of  accurate  theo¬ 
retical  prediction  methods,  a  concurrent  experimental  transonic  research  pro¬ 
gram  on  a  typical  modern  projectile  shape  (Figure  13)  was  carried  out  by  BRL. 
The  primary  experimental  data  consisted  of  surface  pressure  and  turbulent 
boundary  layer  profile  surveys  along  the  aft  portion  of  a  non-spinning  model 
obtained  during  two  different  types  of  tests.  In  the  first  series,  a  pres¬ 
sure-tapped  three  foot  long  sting-mounted  model  was  tested  by  Danberg3'  in 

the  NASA  Langley  Research  Center's  8  foot  transonic  tunnel  at  a  Reynolds 
number  of  13  x  106/M  over  a  Mach  number  range  .8  <  M^  <  .97.  The  boundary 
layer  surveys  were  made  along  both  the  boattail  region  and  some  upstream 
stations  on  the  cylindrical  portion  using  the  sting-mounted  total  head  probe 

shown  in  Figure  14.  In  the  second  test  series38,  a  one  foot  long  model  was 

surveyed  in  the  Naval  Surface  Weapons  Center  White  Oak  Supersonic  Wind  Tunnel 

#2  at  Re^  =  4.5  x  lo6  and  M^  =  .908  using  a  Laser-Doppler  Velocimeter. 

B.  Thin-Layer  Navier-Stokes  Code. 

The  Navier-Stokes  code  with  which  we  have  compared  in  this  study  is  the 
n-invariant  or  Generalized  Axisymmetric  version35.  This  code  solves  the  thin- 
layer  Navier-Stokes  equations  which  are  cast  in  strong  conservation  form.  The 
formulation  allows  for  arbitrary  body  geometries  and  is  solved  using  an 

implicit  approximate  factorization  finite  difference  scheme.  The  "thin-layer" 

approximation  used  requires  Re  »  1;  all  the  viscous  terms  in  the  coordinate 
direction  along  the  body  surface  are  neglected  while  terms  in  the  near  normal 


37.  Danberg ,  J.E.,  Reklis,  R.P. ,  and  Inger,  G.R.,  'Pressure  Distributions  and 
Boundary  Layer  Profiles  on  a  Yawed  Projectile  at  Transonic  Speeds  , " 
University  of  Delaware  Technical  Report  No.  226 ,  Department  of  Mechanical 
and  Aerospace  Engineering,  April  1979  (see  also  AIAA  Paper  79-1551 ,  July 
1979). 

38.  Nietubiaz,  C.J.,  unpublished  wind  tunnel  data.  To  be  published  as  a  BRL 
Memorandum  Report. 
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direction  to  the  body  are  retained.  This  approximation  is  used  because,  due 
to  computer  speed  and  storage  limitations,  fine  grid  spacing  can  only  be  pro¬ 
vided  in  one  coordinate  direction  (usually  taken  as  the  near  normal  direction) 
and  the  grid  spacing  available  in  the  other  two  directions  is  usually  too 
coarse  to  resolve  the  viscous  terms.  The  numerical  algorithm  used  is  a  fully 
implicit,  approximately  factored  finite  difference  scheme  (more  details  of  the 
numerical  method,  algorithm  and  boundary  conditions  can  be  found  in  Reference 
35).  The  numerical  grid  used  for  all  the  Navier-Stokes  computations  is  shown 
in  Figure  15.  The  computational  region  has  been  extended  to  four  model 
lengths  in  front  of  and  four  model  lengths  behind  the  projectile.  The  far 
field  boundary  has  been  set  at  five  model  lengths.  Such  an  extensive  domain 
is  used  to  eliminate  the  possibility  of  any  wave  reflection  back  on  to  the 
model.  The  dark  band  of  lines  near  the  model  surface  results  from  a  cluster¬ 
ing  of  grid  lines  which  are  required  in  order  to  adequately  resolve  the  bound¬ 
ary  layer  region.  The  minimum  spacing  at  the  wall  was  .0002  model  diameters, 
which  resulted  in  at  least  3-5  grid  points  being  in  the  laminar  sublayer.  The 
total  number  of  points  in  the  normal  direction  was  40.  There  were  78  grid 
points  in  the  longitudinal  direction  with  clustering  taking  place  at  X/D  =  3.2 
and  5.3,  near  the  ogive  and  boattail  junction,  respectively. 

C.  Comparisons  of  Theory  and  Experiment 

1.  Pressure  Distributions.  A  typical  comparison  at  =  .95  of  the 

present  composi te  theory  predi cti on  vs.  that  from  the  parabolized  Navier- 
Stokes  code  is  shown  in  Figure  16  to  illustrate  the  acceptable  accuracy  of  the 
inviscid  transonic  small  disturbance  theory  embodied  in  the  former.  Further 
comparisons  between  the  two  methods  and  with  the  experimental  data  are  shown 
in  Figures  17  for  several  other  Mach  numbers,  where  is  it  seen  that  overall 
both  predictions  agree  well  with  the  data.  The  two  rapid  expansion  regions 
followed  by  recompression  shocks,  associated  with  the  body  geometry  disconti¬ 
nuities,  are  clearly  evident  in  the  theoretical  curves  and  are  fairly  well 

modelled  by  both  theories  in  view  of  the  rapid  streamwise  changes  involved. 

The  slight  differences  between  the  two  depend  somewhat  on  the  mesh  sizes  used 
and  also  on  the  fact  that  the  shock  position  of  the  composite  solution  has  not 
been  iterated  with  the  global  boundary  layer  displacement  effect  (as  explained 
above)  whereas  the  Navier-Stokes  solution  automatically  includes  this  effect. 

2.  Boundary  Layer  Velocity  Profiles.  Figure  18  presents  comparisons  for 

a  typical  M^  =  .94  flight  case  of  the  theoretical  and  measured  boundary  layer 
velocity  profiles  at  three  stations  behind  both  the  ogive-cylinder  and  the  aft 
cylinder-boattail  junctions.  Overall,  it  is  seen  that  the  composite  and  thin- 
layer  Navier-Stokes  solutions  are  in  reasonably  good  agreement  with  each  other 
and  the  data  as  regards  the  general  shape  and  streamwi  se-evolution  of  the 

boundary  layer,  although  there  are  some  minor  detailed  differences  which, 

viewed  in  the  light  of  the  customary  experimental  uncertainties  involved, 
appear  minor.  As  expected,  the  greatest  discrepancies  are  found  at  the  most 
aft  station  near  the  boattail  base  where  there  is  very  likely  some  theoretical 
error  due  to  the  crude  base-effect  model  plus  some  possible  experimental  inac¬ 
curacy  due  to  probe  interference  effects  in  the  boattail  interaction  region37. 
Figure  19  presents  analogous  comparisons  illustrating  the  same  conclusions  at 
the  slightly  higher  flight  Mach  number  MOT  =  .97. 
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3.  Displacement  Thickness  Distributions.  Figure  20a  illustrates  the  2-D 
displacement  thickness  growth  along  the  projectile  at  =  .94  predicted  by 
the  composite  theory;  it  indicates  a  significant  interactive  thickening  effect 
across  each  shock  zone,  although  this  is  somewhat  mitigated  (especially  on  the 
boattail)  by  the  thinning  effect  of  the  preceding  sharp  expansion  regions. 
In  Figure  20b  we  compare  with  the  thin-layer  Navier-Stokes  .code  prediction  and 
with  experimental  values  in  the  boattail  region  calculated  from  the  velocity 
profile  data.  The  two  theories  are  in  reasonably  good  agreement  with  experi¬ 
ment  bearing  in  mind  the  joint  analytical  -  experimental  uncertainties  in  the 
aft  part  of  the  boattail.  The  most  notable  theoretical  difference  is  that 
while  the  composite  solution  shows  a  rapid  rise  in  displacement  thickness 
across  the  first  shock  location  and  then  a  gradual  increase  until  the  expan¬ 
sion  of  the  second  corner  is  felt,  the  Navier-Stokes  results  show  a  continuous 
increase  in  displacement  thickness  over  the  cylinder  portion.  This  difference 
is  attributable  to  the  fact  that  the  displacement  thickness  is  difficult  to 
define  in  the  Navier-Stokes  calculation  because  it  does  not  distinguish  the 
boundary  layer  per  se  as  does  the  composite  method. 

Further  comparisons  for  the  lower  flight  Mach  number  case  =  .908  are 

given  in  Figure  21,  showing  similar  features  and  good  experimental  corrobora¬ 
tion  (especially  of  the  composite  theory's  predicted  interactive  6*-  rise)  on 
the  aft  part  of  the  cylindrical  section  as  well  as  on  the  boattail. 

4.  Skin  Friction  Distributions.  Although  no  experimental  values  are 

available  for  comparison,  it  is  nevertheless  instructive  to  examine  the  theo¬ 
retical  predictions  for  this  important  and  sensitive  property  as  a  means  of 

assessing  code  performance  and  the  effects  of  shock- boundary  layer  interaction 
including  regions  of  likely  incipient  shock-induced  local  boundary  layer  sepa¬ 
ration.  Indeed,  the  prediction  of  skin  friction  behavior  is  usually  the  most- 
demanding  test  of  any  theoretical  method;  this  is  particularly  true  in  the 
present  application  because  of  the  two  regions  of  rapid  expansion/compression 
along  the  body.  Typical  results  are  illustrated  in  Figure  22  for  the  Mro  = 

.908  case,  where  it  can  be  seen  that  both  the  composite  and  Navier-Stokes 
solutions  indicate  a  sharp  Cf  rise  in  the  expansion  zones  at  both  corners 

followed  immediately  by  rapid  Cf  decrease  in  the  shock-interaction  zone  with  a 

subsequent  downstream  recovery.  The  composite  solution  predicts  noticeably 
larger  peak  values  of  these  maxima  and  minima  then  does  the  Navier-Stokes 
code;  however,  part  of  this  difference  is  seen  to  be  due  to  the  latter  s' 

underprediction  of  the  upstream  turbulent  skin  friction  level  along  the  ogive 
section;  when  this  discrepancy  is  corrected  for,  the  two  methods  are  seen  to 
be  in  good  quantitative  (as  well  as  qualitative)  agreement,  especially  consi¬ 
dering  the  severe  Cf  changes  involved. 

Similar  conclusions  may  be  drawn  from  the  predictions  illustrated  in 
Figures  23  and  24  for  =  .94  and  .97,  respectively.  In  all  these  results  it 
can  be  observed  that  the  shock-boundary  layer  interaction  zones  noticeably 
reduce  Cf  well  downstream  of  the  shock.  Moreover,  this  effect  is  especially 

strong  in  the  boattail  region  where  shock-induced  local  separation  appears 
immi nent. 
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D.  Predicted  Interaction  Effects  in  the  Presence  of  Spin 


Although  experimental  data  and  Navier-Stokes  calculations  for  spinning 
bodies  at  transonic  speeds  are  presently  not  available,  it  is  nevertheless 
interesting  to  examine  the  influence  of  spin  on  the  interactive  flow  field 
predicted  by  the  composite  viscous-inviscid  solution.  Accordingly,  the 
effects  of  increasing  spin  rate  ft  were  studied  for  a  typical  full  scale  pro¬ 
jectile  flight  case;  the  results  are  presented  in  Figures  25a-c  for  the 
streamwise  distributions  of  boundary  layer  shape  factor,  displacement  thick¬ 
ness  and  skin  friction  as  a  function  of  ft.  It  is  seen  that  spin  rates  of  10 
to  15  thousand  RPM  or  higher  exert  a  modest  but  discernible  influence  that 
enhances  the  interactive  thickening  on  both  the  cylinder  and  boattail  regions, 
along  with  an  attendant  reduction  in  shape  factor  and  some  increase  in  Cf  in 

the  upstream  ogive  region.  It  appears  questionable,  however,  whether  the 
magnitude  of  these  spin  effects  on  the  interaction  could  be  distinguished 
within  the  typical  uncertainty  level  of  experimental  measurement. 


VI.  CONCLUDING  REMARKS 

This  report  has  presented  a  detailed  triple-deck  theory  of  the  transonic 
shock-nonseparating  turbulent  boundary  layer  interaction  zones  on  spinning 
axisymmetric  bodies  at  zero  angle  of  attack.  The  theory  allows  a  rather 
general  non-equilibrium  turbulent  boundary  layer  history  upstream  of  the 
interaction.  We  further  described  the  application  of  this  theory  as  a  local 
interactive  model  imbedded  in  a  global  composite  viscous-inviscid  transonic 
flow  field  prediction  method  which  can  handle  two  or  more  such  interaction 
zones.  Detailed  comparisons  with  both  experimental  data  and  thin-layer 
Navier-Stokes  code  predictions  were  presented  to  verify  the  accuracy  and 
reliability  of  the  composite  method  solutions  for  pressure,  displacement 
thickness,  velocity  profiles,  shape  factor  and  skin  friction  distributions 
along  the  body  over  a  wide  range  of  Mach/Reynolds  number  conditions. 

Desirable  improvements  in  the  present  work  would  include  extension  to  the 
case  of  small  (a  <  3°-4°)  angles  of  attack  and  better  treatment  of  the  aft 
region  near  the  base  to  account  for  the  upstream  influence  and  true  flow 
conditions  there. 
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Figure  1.  Typical  Double-Shock  Flow  Pattern  on  a  Transonic  Projectile 
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Figure  2.  Separated  versus  Non-Separati  ng  Shock-Turbulent 
Boundary  Layer  Interaction  Patterns 
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Figure  4.  Effective  Wall  Concept  Associated  with 
the  Inner  Deck  Solution  (Schematic) 
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Figure  5.  The  Displacement  Thickness  Function  H(t) 
and  Disturbance  Shear  Function  S(t) 
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a.  Typical  Pressure  Data 
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Figure  6. 


Comparison  of  Tripie-Deck  Theory  with  the  Experimental 
Data  of  Ackeret,  Feldmann  and  Rott 
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Figure  7.  Comparison  of  the  Predicted  Interactive  Shape 
Factor  Variation  with  ONERA  Data 
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Figure  8.  Comparison  of  Triple-Deck  Theory  with  DFVLR-AVA 
(Go)  Measurements  on  a  Supercritical  Airfoil 


46 


Figure  10.  Schematic  of  the  Expansion-Corner  Regions 
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Figure  11.  Shock  Obliquity  Due  to  Viscous  Effects 
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Figure  13.  Typical  Projectile  Model  for  Experimental  Studies 


Figure  14. 


Sting-Mounted  Total  Head  Probe 


in  Danberg's  Experiment 


50 


— 


-  I  I  i  i  i  i  |  i  i  i  i  j'  i  i  i  i  j  r  r  r  i  |  r  i  i  i  |  I  i  I  I  | 

0  L  2  3  *  5  6 

X/D 

Figure  15.  Grid  System  Used  in  the  Navier-Stokes 

Calculations  of  Nietubicz  et  al  (Ref.  35) 
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Figure  16.  Comparison  of  Transonic  Small  Disturbance  Theory  and  Navier- 
Stokes  Predictions  of  the  Axial  Pressure  Disbribution 
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Figure  17,  Comparison  of  Composite  Invi scid-Vi scous  and  Navier-Stokes 
Pressure  Distributions  versus  Experiment 
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Figure  18.  Comparison  of  Composite  Solution,  Navier-Stokes  and 
Experimental  Boundary  Layer  Profiles  at  Mm  =  .94 
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Figure  19.  Comparison  of  Composite,  Navier-Stokes  and  Experimental  Boundary  Layer  Profiles  at  M 
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b.  Comparison  with  Navier-Stokes  Results  and  Experimental  Data 


Figure  20.  Displacement  Thickness  Distributions  at  =  .94 
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Figure  22.  Composite  Interaction  Theory  and  Navier-Stokes  Predictions  for  Skin  Friction 
Distribution  at  M  =  .908 
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Figure  24.  Theoretical  Skin  Friction  Predictions  at  =  .97 
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Figure  25.  The  Typical  Influence  of  Spin  on  the  Viscous-Inviscid 
Interaction  Field  Predicted  by  the  Composite  Theory 
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LIST  OF  SYMBOLS 


A 

Van  Driest-Cebeci  wall  turbulence  damping  parameter 

cf 

skin  friction  coefficient,  2xw/pe  U0  2 

eo  0 

CP 

pressure  coefficient,  2  p'/Pe  2 

0  0 

H 

boundary  layer  shape  factor,  S*/0* 

Hi 

incompressible  shape  factor 

M 

Mach  number 

P 

static  pressure 

P' 

interactive  pressure  perturbation,  p-p-| 

Ap 

pressure  jump  across  incident  shock 

Rer  Re6 

Reynolds  number  based  on  length  i  and  boundary  layer  thickness, 
respecti vely 

T 

absolute  temperature 

T 

basic  interactive  wall-turbulence  parameter  (see  Eq.  23) 

u',  v1 

streamwise  and  normal  interactive  disturbance  velocity  components, 
respectively 

u0 

undisturbed  incoming  boundary  layer  velocity  in  x-di recti  on 

X,  y 

streamwise  and  normal  distance  coordinates  (origin  at  the  inviscid 
shock  intersection  with  the  wall) 

Yweff 

effective  wall  shift  seen  by  interactive  inviscid  flow  (see  Fig.  3) 

3 

>/Mi2-T 

Y 

specific  heat  ratio 

6 

boundary  layer  thickness 

6* 

boundary  layer  displacement  thickness 

6SL 

inner  deck  sublayer  thickness 

eT 

kinematic  turbulent  eddy  viscosity 
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LIST  OF  SYMBOLS  (continued) 


£j‘  interactive  perturbation  of  turbulent  eddy  viscosity 

i  y/«0 

y  molecular  coefficient  of  viscosity 

v  y/p 

“  viscosity-temperature  dependence  exponent,  y  ~  Tw 

p  densi ty 

0*  boundary  layer  momentum  thickness 

t  total  shear  stress 

t1  interactive  perturbation  of  total  shear  stress 

Subscripts 

1 
e 

inc 
inv 
o 


undisturbed  inviscid  values  ahead  of  incident  shock 
conditions  at  the  boundary  layer  edge 
incompressible  value 
inviscid  disturbance  solution  value 
undisturbed  incoming  boundary  layer  properties 
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COMPOSITE  LAW  OF  THE  WALL- LAW  OF  THE  WAKE 
TURBULENT  VELOCITY  PROFILE  RELATIONSHIPS 
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APPENDIX  A 


COMPOSITE  LAW  OF  THE  WALL-LAW  OF  THE  WAKE 
TURBULENT  VELOCITY  PROFILE  RELATIONSHIPS 


Because  of  its  convenient  analytical  form,  accurately  blended  representa¬ 
tion  of  the  combined  Law  of  the  Wall-Law  of  the  Wake  behavior  and  generality, 
we  have  adopted  Walz's  model32  for  the  incoming  turbulent  boundary  layer  up¬ 
stream  of  the  interaction.  For  the  low  Mach  number  adiabatic  wall  conditions 
appropriate  to  transonic  applications,  it  may  be  satisfactorily  corrected  for 
compressibility  effects  by  the  Eckert  Reference  Temperature  method  (which 
under  these  conditions  is,  in  fact,  comparable  in  accuracy  to,  but  far  simpler 
to  implement  than,  the  Van  Driest  compressibility  transformation  approach). 

Let  tt  be  Coles  (incompressible)  Wake  Function,  n  =  y/6  and  denote  for 

convenience  R  =  .41  Re6Q*/[(l  +  "HTw/Te)1  +  ^  where  Tw/Te  “  1  +  -18  Me?  and 

to  “  .76  for  a  perfect  gas;  then  the  compressible  form  of  Walz's  composite 
profile  may  be  written 

■jj“  =  1  +  — y/— 2“^Te  ^  +  ^  11 2  (8~2n)  -  2ir 

oe 

( A— 1 ) 

+  in  (^2)  -  (2.15  +  1.235  Rn)e_,3Rn] 


subject  to  the  following  condition  linking  it  to  Cfo  and  Re6Q*  that  derives 
from  the  u  ue  matching  condition  at  the  boundary  layer  edge: 


2 it  +  2.15  +  In  (1+R)  = 


.41 


(A-2) 


Eqs.  (A-l)  and  (A-2)  have  the  following  desireable  properties:  (a)  for 
n  >  .10  or  so,  Uo  is  dominated  by  a  Law  of  the  Wake  behavior  which  correctly 
satisfies  both  the  outer  limit  conditions  Uo/Ue  1  and  dU  /dy  0  as  n  1; 

(b)  on  the  other  hand,  for  very  small  n  values,  Uo  assumes  a  Law  of  the  Wall- 
type  behavior  consisting  of  a  logarithmic  term  that  is  exponentially  damped 
out  extremely  close  to  the  wall  into  a  linear  laminar  sublayer  profile  U0/Ue 

-  Rn  as  n  0;  (c)  Eq.  (A-l)  may  be  differentiated  w.r.t.  n  to  yield  an 
analytical  expression  for  dUQ/dy,  which  proves  advantageous  in  solving  the 
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middle  and  inner  deck  interaction  problems  (see  text)  where  dMQ/dy  must  be 
known  and  vanish  at  the  boundary  layer  edge. 


The  use  of  the  incompressible  form  of  (A-l)  in  the  defining  integral 
relations  for  6..*  and  0>  yields  the  following  relationship  that  links  the 

wake  parameter  to  the  resulting  incompressible  shape  factor  =  ( 6^ */0^ *) : 


1+1 ,59tt  +  .75  tt  2 

1  +  7T 


(A-3) 


Eqs.  (A-2)  and  (A-3)  together  with  the  defining  relation  for  R  enable  a 
general  and  convenient  parameterization  of  the  profile  (and  hence  the  interac¬ 
tion  that  depends  on  it)  in  terms  of  the  shock  strength  Me^  the  local  dis¬ 
placement  thickness  Reynolds  number  Re^*,  and  one  additional  physical  para¬ 
meter  upstream  of  the  shock.  This  parameter  is  taken  to  be  the  shape  factor 
H.jj  because  it  reflects  the  upstream  history  of  the  incoming  boundary  layer 

including  possible  pressure  gradient  and  surface  mass  transfer  effects.  With 
this  prescribed,  the  aforementioned  three  equations  may  be  solved  simulta¬ 
neously  for  the  attendant  skin  friction  value  Cf0,  the  value  of  R  and,  if 

desired,  the  n  value  appropriate  to  these  flow  conditions. 
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